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FOREWORD 


The  work  performed  under  Contract  DAAK40-78-C-0031  that  resulted  in 
specific  deliverables  (principally  computer  programs)  is  described  in  two 
volumes  of  a  final  report.  This  is  Volume  2,  which  describes  the  work  that 
led  to  the  extended  target  model.  Volume  1,  entitled  "Distributed  Clutter 
Model:  Real-Time  FFT  Program,"  describes  the  work  relevant  to  that  subject. 

The  AP120B  array  processor  programs  described  in  this  report  (but  docu¬ 
mented  elsewhere)  were  designed  and  developed  by  Dr.  I.  P.  Bottlik,  as  well 
as  the  Interfaces  to  the  host  computer  and  Digital  Signal  Generator. 

Dr.  Bottlik  was  assisted  by  Mr.  D.  L.  Brandon. 
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1.  INTRODUCTION  AND  SUMMARY 


This  report  describes  a  real-time  simulation  procedure  for  generating 
replicas  of  an  extended  target  signal  to  be  radiated  from  the  RFSS  array 
to  a  missile  system  radar  under  test.  The  procedure  is  purely  deterministic. 
It  is  based  on  identifying  in  the  model  the  location  and  amplitude  of  the 
predominant  scattering  centers  on  the  target;  during  the  real-time  simulation, 
rays  are  traced  to  each  scatterer  so  that  the  radiated  signal  is  the  phasor 
superposition  of  the  component  rays.  The  missile  and  target  geometries  are 
modeled  with  six  degrees  of  freedom,  and  the  target  motions  of  roll,  pitch, 
and  yaw  are  included.  The  instantaneous  target  aspect  is  determined  so  that 
any  scattering  dependencies  on  target  aspect  can  be  Included  in  the  model. 

The  range  dimension  of  the  target  is  sorted  into  range  bins  corresponding 
to  taps  on  a  tapped-delay  line.  For  each  tap  an  angular  glint  centroid  is 
computed.  - 

-  v  The  procedure  permits  any  target  to  be  simulated;  the  scattering  centers 
on  the  target  must  be  identified  beforehand  as  part  of  the  model.  This  in¬ 
formation  forms  the  data  base  from  which  the  signal  is  generated.  Listings 
of  a  Fortran  program  are  included  in  Appendix  B.  As  written,  the  program 
will  not  run  in  real  time.  It  is  assumed  that  some  portion  of  it  will  be 
Implemented  on  an  AP120B  array  processor  to  satisfy  the  real-time  require¬ 
ment. 
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2.  THE  MEDIUM-RANGE  TARGET  MODEL 

Three  types  of  extended  target  models  were  discussed  in  Reference  1, 
the  short-*  medium-,  and  long-range  models.  The  long-range  model  has  only 
limited  application  in  the  RFSS,  and  the  short-range  model,  which  is  based 
on  radiating  signals  into  specific  angles  on  the  receive  beam,  requires  ac¬ 
curate  knowledge  of  the  receive  antenna  boresight  axis  position.  Because 
of  these  disadvantages  it  was  decided  to  concentrate  on  the  medium-range 
model,  which  assumes  that  one  is  always  working  in  the  linear  portion  of 
the  monopulse  receive  beam.  The  target  model  is  uncoupled  from  the  receive 
antenna  position  as  long  as  this  assumption  is  valid.  In  this  section  we 
derive  the  form  of  the  modulation  signals  for  the  medium-range  model. 

2.1  The  Scintillating  Target 

For  the  moment,  let  us  assume  that  the  target  is  composed  of  a  number 
of  point  scatterers  at  the  same  effective  range  (no  range  extent  compared 
to  the  resolution  cell,  but  the  scatterers  may  differ  in  range  by  many  wave¬ 
lengths).  If  denotes  the  radar  cross  section  (RCS)  of  the  kcn  scatterer 

as  viewed  from  a  particular  target  aspect,  then  we  can  assign  a  complex 
cli  2 

quantity  to  the  ktn  scatterer,  V^,  where  |V^|  -  and  the  phase  accounts 

for  the  slightly  different  range  (measured  in  wavelengths)  among  the  scat¬ 
terers  (or  differences  in  the  reflection  properties) .  If  all  scatterers 
are  illuminated  with  the  same  transmit  and  receive  gain  we  can  replace  the 
set  of  scatterers  by  a  single  scatterer  with  the  effective  complex  reflec¬ 
tion  coefficient  given  by 


(1) 


The  phasor  summation  in  (1)  accounts  for  the  target  signal  scintillation  as 
measured  by  the  radar.  With  the  exception  of  a  scale  factor  to  account  for 
the  radar  range  equation,  it  is  the  modulation  function  that  is  applied  to 
the  delayed  transmitted  waveform  to  simulate  the  target. 
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2.2  The  Glint  Centroid 


Let  us  take  the  same  target  model  above  and  assume  that  the  receiver 
is  a  monopulse  system  so  that  it  is  capable  of  measuring  an  angle  to  the 
target.  We  will  again  assume  that  the  target  has  no  range  extent.  If  all 
scatterers  are  illuminated  by  the  same  transmit  antenna  gain,  and  all  are 
vithin  the  linear  region  of  the  receive  monopulse  beam,  the  radar  will 
measure  an  angle  to  the  target  that  is  effectively  the  electrical  centroid 
given  for  the  one-dimensional  case  by 


0  +  Re 
o 


S<V9o>Vk 


(2) 


where  0q  is  a  reference  angle  and  6^  is  the  angle  associated  with  the 
scatterer. 

Equation  (2)  is  insensitive  to  the  angle  of  the  receive  antenna  bore- 
sight  axis.  It  is  valid  as  long  as  all  scatterers  are  within  the  linear 
region  of  the  monopulse  beam  (i.e. ,  the  target  model  will  be  uncoupled  from 
the  radar  antenna).  In  practice,  this  linear  region  extends  approximately 
only  to  about  the  half-power  beamwidth  of  the  sum-channel  beam.  In 
Appendix  A  we  derive  the  maximum  angular  extent  of  the  target  for  which 
the  above  assumptions  are  valid,  and  it  is  given  roughly  by  the  half-power 
beamwidth  (one-way) . 


2.3  Range  Extent 

Modern  radars  usually  have  high  range  resolution  so  that  the  target  may 
extend  over  several  range  resolution  cells  (this  is  the  definition  of  a 
range-extended  target).  In  this  case  only  a  few  of  the  scatterers  on  a  tar¬ 
get  will  fall  within  any  given  range  cell.  The  response  in  a  range  cell  is 
given  by  a  slightly  modified  from  of  (1)  and  (2)  as 
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0-0  +  Re 
o 


i-Wic 
k _ 

K\ 


(4) 


where  accounts  for  the  variable  range  gate  weighting  of  the  ktil  scatterer 
in  the  range  cell  of  interest.  It  is  a  function  of  only  the  differential  de¬ 
lay  between  the  k^  scatterer  and  the  center  of  the  range  gate  (see  Section  4 ) . 

In  general,  a  strict  application  of  (3)  and  (4)  requires  that  the  range 
gate  timing  be  known  in  the  receiver.  However,  as  was  derived  in  Reference  1 
it  is  possible  to  accurately  simulate  the  range  extent  of  the  target  when 
the  timing  is  unknown  by  the  use  of  a  resampling  technique  based  upon  a  tapped- 
delay  line  implementation.  The  taps  must  be  spaced  no  further  apart  than 
50Z  of  the  range  resolution  cell  size.  The  technique  consists  simply  of  ap¬ 
plying  (3)  and  (4)  to  each  tap  of  the  tapped-delay  line.  The  weights  are 
thus  defined  as  the  tap  weights,  and  they  are  a  function  of  the  differential 
delay  between  the  location  of  the  scatterer  and  the  tap.  For  a  spacing  of 
50Z  of  the  range  resolution  cell  (the  most  efficient  spacing)  only  four  taps 
will  receive  any  significant  contribution  from  each  scatterer  (see  Section  4). 

2.4  Implementing  the  Medium-Range  Model  on  the  RFSS 

The  use  of  (3)  and  (4)  creates  a  sequence  of  effective  complex  voltage 

A 

samples  V^(n)  and  glint  centroids  0(n)  as  a  function  of  range  as  is  sketched 
in  Figure  1.  In  order  to  simulate  the  composite  signal  we  have  to  create  a 


AR  -  Tap  Spacing 


Figure  1.  Physical  Arrangement  of  Resampled  Scatterers  for  Range-Extended 
Target 
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different  glint  angle  for  each  tap  on  a  tapped -delay  line.  One  approach 
would  be  to  utilize  a  separate  RFSS  signal  generation  channel  for  each  tap, 
but  the  approach  would  be  limited  to  four  taps  because  there  are  only  four  ' 

RFSS  channels.  A  far  better  approach  would  be  to  utilize  the  supertriad  de¬ 
vised  by  Mott  [2],  which  we  now  discuss. 

To  simulate  a  point  target  in  the  RFSS,  signals  are  directed  to  three 
horns  on  the  array  (a  triad)  and  with  proper  amplitude  and  phase  control  on  * 

each  horn  the  point  target  can  be  made  to  appear  anywhere  within  the  triad. 

For  our  purposes  we  will  not  restrict  the  three  sources  to  coincide  with 
the  ABC  triad  configuration  of  the  RFSS  array.  We  will  use  separate  RFSS 
channels  to  generate  each  signal  so  that  the  three  points  of  radiation  can 
be  placed  anywhere  on  the  RFSS  array.  We  will  designate  this  configuration 
as  a  "supertriad."  The  elements  of  the  supertriad  can  be  spaced  further 
apart  than  the  horn  spacing  in  the  conventional  ABC  triad. 

For  a  range-extended  target  we  will  build  upon  the  supertriad  with  a 
tapped-delay  line  as  shown  in  Figure  2.  The  RFSS  signal  channels  are  desig¬ 
nated  A,  B,  and  C  in  the  figure  to  coincide  with  the  elements  of  the  super¬ 
triad.  The  A,  B,  and  C  channel  signals  will  each  be  radiated  from  fixed  ^ 

points  on  the  array,  regardless  of  the  delay  (until  the  engagement  geometry 

is  updated) ;  however ,  by  controlling  the  modulation  signals  (A1 ,  A2 . 

AS,  Bl,  B2,  ....  C8  in  Figure  2)  the  phase  center  can  be  made  to  appear 

anywhere  within  the  supertriad  (and  slightly  outside  it)  at  each  tap  on  the  ( 

tapped-delay  line.  Thus  we  have  the  situation  that  is  sketched  in  Figure  3. 

Note  that  the  location  of  the  phase  center  in  angle  at  each  tap  can  be  chosen 
Independently  of  the  other  taps. 

2.5  The  Modulation  Signals  at  the  Taps  ( 


Equations  (3)  and  (4)  define  the  signal  amplitude,  phase,  and  angle 
(there  will  be  two  angles  involved)  corresponding  to  one  tap  of  the  tapped- 
delay  line.  To  compute  the  modulation  signal  that  is  applied  to  the  super¬ 
triad,  let  us  refer  to  Figure  4  where  we  place  the  origin  of  the  x,y  coordi¬ 
nate  system  at  the  centroid  of  the  supertriad  composed  of  equal  sides.  In 
other  words.  Point  A  is  located  at  the  coordinate  (0,  d//3).  Point  B  at 


INPUT  SIGNAL 


D  -  DELAY 
M  -  MODULATOR 
£  «  SUMMER 


CHANNEL  A  CHANNEL  B 


Figure  2.  Supertriad  Configuration  with  Delay  Line 
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The  reference  angles  (8  ,  4>  )  coincide  with  the  origin  of  the  x,y  coordi- 

o  o 

nates.  The  modulation  signals  that  are  applied  to  the  elements  of  the 
supertriad  are  now 

V.  -  V  4  +  ~  f)  (7) 

A  e  3  ^  d 


V  (4  +  f  -  —  J) 
6  3  d  /J  d 


VC  Ve  4  d  ^  d> 


The  above  computations  must  be  applied  to  each  tap. 

2.6  Choosing  the  Size  of  the  Supertriad 

In  general,  the  size  of  the  supertriad  should  be  chosen  to  encompass 
the  angular  extent  of  the  target.  At  long  range  we  can  choose  d  small  so 
the  accuracy  will  be  best  (but  d  should  not  be  less  than  the  RFSS  horn 
spacing).  As  the  target  range  decreases  we  must  increase  d  accordingly. 

As  long  as  the  elements  of  the  supertriad  are  within  the  linear  region  of 
the  monopulse  beam,  there  will  be  no  error  due  to  the  fact  that  we  are 
radiating  from  three  points  simultaneously,  instead  of  a  single  one  at  the 
desired  phase  center.  However,  we  will  eventually  reach  a  limiting  value 
of  d  where  we  are  no  longer  within  the  linear  region  of  the  monopulse  beam. 

To  determine  this  limit  let  us  approximate  the  one-dimensional  response 
to  the  monopulse  system  by 
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«b«re  8  is  the  angle  of  a  point  acatterer  measured  from  the  boresight  axis. 
Let  us  place  two  horns  spaced  by  A  -  d/D  as  in  Figure  5  with  one  at  an  angle 
®o  ^Tom  *ke  origin,  and  we  will  place  a  point  scatterer  at  the  origin.  From 
the  left  horn  we  will  radiate  a  signal  with  a  relative  amplitude  8/A  and 
from  the  right  horn  a  signal  of  relative  amplitude  (1  -  8/A).  The  composite 
response  of  the  monopulse  system  is  the  sum  of  the  individual  responses  as 

6  -  8(l-a82)(l-8/A)  +  © -A) [l-a(8-A)2]  (8/A) 


-  a6(A-8)(A-26)  (11) 

The  angle  would  be  zero  if  the  system  were  linear  (i.e. ,  a  ■  0) ,  so  (11) 
represents  the  error  in  simulating  a  scatterer  on  the  boresight  with 
two  horns  configured  as  in  Figure  5.  A  local  maximum  (extremum)  of  (11) 
occurs  at  8  •  A (1/2  +  1//12)  and  is  given  by 


ext 


3^12 


ctA“ 


(12) 


A  typical  value  for  a  is  1.70/82^,  where  8^^  is  the  one-way  half¬ 
power  beamwidth  of  the  monopulse  sum  channel.  For  this  case  and  8.  »  12° 

3dB 

we  have  plotted  the  error  given  by  (11)  in  Figure  6.  The  abscissa  is 
®  ”  A/2,  so  that  a  value  of  8  —  A/2  m  0  corresponds  to  the  two  horns  being 
placed  symmetrically  about  the  origin.  The  curves  are  also  symmetrical 
about  the  origin  in  Figure  6.  For  a  -  1.70/8^  we  can  rewrite  (12)  as 


(13) 
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Figure  5.  Simulating  Point  Scatterer  with  Two  Horns 


3.  ENGAGEMENT  GEOMETRY  AND  COORDINATE  TRANSFORMS 

In  this  section  we  define  the  quantities  that  are  required  to  specify 
the  engagement  geometry  and  the  transformations  used  to  obtain  the  ray  co¬ 
ordinates  of  each  scatterer  in  the  missile  coordinate  system. 

3.1  Specification  of  Target  Position  and  Dynamics 

At  a  given  instant  of  time  the  target  c.g.  is  characterized  by  a  loca¬ 
tion  (x  ,  y  ,  z  )  and  rate  (&  ,  f  ,  i  )  in  an  inertial  coordinate  system 
o  o  o  o  o  o 

referenced  to  the  ground.  The  x-y  plane  is  parallel  to  the  ground  and  the 

z-axis  is  positive  downward. 

The  target  coordinate  system  is  sketched  in  Figure  7  and  is  defined  such 
that  the  x-axis  is  the  longitudinal  axis  of  the  aircraft  target,  positive  in 

the  direction  of  the  nose;  the  y-axis  is  in  the  direction  of  the  right  wing; 

and  the  z-axis  is  down.  The  origin  of  the  target  coordinate  system  is  at 
the  target  c.g. 

The  orientation  of  the  target  is  given  by  a  series  of  three  consecutive 
coordinate  system  rotations,  the  order  of  which  is  Important.  First,  the 
target  is  oriented  so  that  its  x,y,z  axes  are  parallel  to  the  ground- 
referenced  x,y,z  axes.  Then  the  following  rotations  (clockwise,  looking 
out  from  the  coordinate  origin)  are  applied: 


Figure  7.  Definition  of  Target  Coordinate  System 
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1.  a  rotation  ¥  about  the  z-axis 

2.  a  rotation  @  about  the  y-axis 

3.  a  rotation  4  about  the  x-axls 

The  target  is  under going  roll,  pitch,  and  yaw  motion  in  its  own  coordi¬ 
nate  system.  We  define  P,Q,R  to  be  the  roll,  pitch,  and  yaw  body  rotation 
rates  (clockwise  looking  out  from  the  coordinate  origin)  as  sketched  in 
Figure  8.  The  relationship  of  P,Q,R  to  the  angles  ¥,0,$  is 


0  ■  Q  cos  4  -  R  sin  4 

(14) 

♦  ■  P  +  (Q  sin  $  +  R  cos  4) 

tan  0 

(15) 

• 

¥  -  (Q  sin  4  +  R  cos  4)  sec 

0 

(16) 

In  general,  the  excursion  of  these  motions  will  be  over  small  angles. 


Figure  8.  Roll,  Pitch,  and  Yaw  Motions 

3.2  Specification  of  Missile  Position  and  Dynamics 

A  the  same  instant  of  time  for  which  the  target  dynamics  are  defined, 

we  define  the  location  of  the  missile  as  (x  ,y_»z  )  and  the  velocity  as 

m  m  m 

(4  ,y  ,i  ),  where  the  coordinate  system  is  the  same  ground-referenced  system 

8  i  m 

/  used  for  the  target. 


3,3  Specification  of  Target  Scatterers 


The  target  aspect  is  defined  in  Figure  9,  where  the  azimuth  and  ele¬ 
vation  angles  of  the  radar  (as  viewed  from  the  target)  are  a  and  g,  respec¬ 
tively.  The  aspect  angle  from  the  x-  or  roll-axis,  which  we  designate  as  y, 
is  given  by 

cos  y  »  cos  a  cos  6  (17) 


Figure  9.  Target  Aspect 


At  a  particular  target  aspect,  a  set  of  scattering  centers  will  be  speci¬ 
fied  in  terms  of  the  location  (x^.y^.z^)  and  amplitude  (A^)  of  each  scat- 
terer,  where  the  subscript  k  designates  the  scatterer  number  and  amplitude 
is  defined  as  the  square  root  of  radar  cross  section  (RCS).  In  general, 
A^,x^,y^,z^  will  all  be  functions  of  the  target  aspect.  A  specification  of 
thia  functional  relationship  is  the  target  model. 


3.4  The  ABC  Vectors 

Let  us  define  three  mutually  perpendicular  unit  vectors,  the  origin  of 
which  is  at  the  target  c.g.  Vector  A  points  to  the  missile  radar.  Vector  B 
is  in  the  horizontal  plane  pointing  to  the  left  as  viewed  from  the  missile, 
and  Vector  C  completes  the  definition  of  a  right-hand  system  (it  points 
downward  as  viewed  from  the  missile) .  In  terms  of  the  conventional  radar 


16 


cootdlnates  (referenced  Co  the  ground) ,  A  points  in  the  direction  of  de¬ 
creasing  range,  B  in  the  negative  azimuth  direction,  and  C  in  the  negative 
elevation  direction.  In  terms  of  the  RTSS  array,  A  points  from  the  target  c.g. 
on  the  RFSS  array  to  the  missile  radar,  B  points  left  as  viewed  from  the 
missile,  and  C  points  down. 

From  the  previous  definitions  we  have  the  direction  cosines  of  A  in  the 
ground  referenced  coordinate  system  given  by 


o 

u 

✓-s 

H° 

1 

z* 

1 

wM 

(18) 

• 

1 

o'4 

«■ 

o 

(19) 

a  »  (z  -  z  )/r 
z  m  o'  o 

(20) 

vhere  r  is  the  range  to 
0 

the  target  c.g.  given  by 

ro  “  (xm 

2  2  2 
-x)  +(y  -  y  )*  +  (z  -  z) 

o  Jo  m  o 

(21)' 

Let  us  define 

p2 

-  (x  -  x  )2  +  (y  -  y  )2 
no  n  o 

(22) 

so  that  we  can  write  the 

direction  cosines  of  the  B  and  C  vectors  as 

\  -  -  <yffl  -  yQ)/p 

(23) 

by  "  (xm  "  Xo)/p 

(24) 

bz-o 

(25) 

H 

1 

1 

0N 

1 

N 

O 

v-/ 

(26) 

V  ■  V1.  -  *.> 

(27) 

cz  -  p 

(28) 

Mote  that  C  »  A  x  B. 

3.5  Transformation  of  a 

Vector  to  Target  Coordinates 

Given  an  arbitrary  vector  A  in  the  ground-referenced  coordinate  system, 
we  can  define  it  in  the  target  coordinate  system  by  three  successive  coordi- 
*  oate  rotations.  First,  if  we  rotate  the  2-axis  by  V  (clockwise  looking  out) 
we  obtain  the  direction  cosines  in  the  new  system 
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u  *  a  cos  4f  +  a„  sin  ¥ 

XX  7 

u  •  -a  sin  ¥  +  a  cos  f 

y  x  7 

u  ■  a 
z  z 

Next,  let  us  rotate  the  7-axis  by  9  (clockwise  looking  out)  so  that 


v  ■  u  cos  0  -  u  sin  9 

XX  z 


u  sin  0  +  u  cos  0 

X  2 


(29) 

(30) 

(31) 

(32) 

(33) 

(34) 


Finally,  we  will  rotate  the  x-axis  by  <5  (clockwise  looking  out)  to  obtain 


w  »  V  cos  $  +  V  sin  $ 
7  7  * 

wz  *  ~v„  ain  «  +  V  cos  4 
Z  7  z 


(35) 

(36) 

(37) 


Since  several  vectors  will  be  transformed  in  this  maimer  we  will  desig¬ 
nate  the  resulting  vector  in  the  target  coordinate  system  as  W(A) ,  where  the 
original  vector  was  A.  The  direction  cosines  of  W(A)  are  designated  as 

W  (A) ,  w  (A),  and  w  (A), 
x  y  z 

The  A.B.  and  C  vectors  (defined  in  Section  3.4)  in  the  target  coordinate 
system  are  thus  W(A),  W(B),  and  ff(C).  Note  that  W(C)  -  W(A)  x  w(B),  so  that 
after  transforming  the  first  two  vectors  we  could  find  the  third  by 

(38) 

(39) 

(40) 


*  (C)  ■  w  (A)w  (B)  -  w  (A)w  (B) 
*  /  Z  Z  / 


w  (C)  -  w  (A)w  (B)  -  w  (A)w  (B) 
y  z  x  x  z 

"(C)  -  w  (A)w  (B)  -  w  (A)w  (B) 
z  *  7  7  x 

3.6  Target  Aspect  Angle 


Let  us  refer  to  Figure  9  where  W(A)  is  the  vector  pointing  to  the  radar. 
The  azimuth  angle  a  is  given  by 


ten  a  •  v  (A)/w  (A) 
7  x 


(41) 


and  Che  elevation  angle  $  by 


sin  3  -  -w2  (A)  (42) 

(remember  chat  W(A)  is  a  unit  vector).  The  aspect  angle  measured  from  the 
roll  axis ,  y ,  is  given  by 

cos  y  ■  cos  a  •  cos  8  ■  w^(A)  (43) 

3.7  Projection  of  Scatterers  onto  ABC  Vectors 

Given  a  scattering  center  located  at  in  target  coordinates, 

the  projection  of  this  coordinate  onto  the  A,B,C-vectors  is  simply  the  dot 
product  onto  each  vector.  Thus  ve  will  write 

Aak  "  Vx®  +  Yy^  +  ZkWz  ^  (44) 

dbk  -  XjW^CB)  +  ykwy (B)  +  zkwz (B)  (45) 

Ack  “  XjWx(C)  +  7kwy(C)  +  zkwz(C)  (46) 

where  we  are  interpreting  the  scatterer  location  to  be  an  increment  applied 

to  the  target  c.g.  While  the  distance  from  the  radar  to  the  scatterer  is 
based  on  the  sum  of  the  vectors  from  the  radar  to  the  target  c.g.  and  the 
target  c.g.  to  the  scatterer,  we  can  considerably  simplify  matters  by  as¬ 
suming  that  the  latter  vector  is  small  in  comparison  with  the  first  (such 
an  assumption  is  consistent  with  the  medium-range  target  model) .  Thus  the 
range  to  the  scatterer  is  essentially  r^-Aa^;  thus  the  differential  range  is 
defined  by 

Arfc  -  (47) 

Viewed  from  the  radar  the  azimuth  and  elevation  angles  (measured  from  the 
target  c.g.)  are  -Abk/rQ  and  -Ack/rQ,  respectively,  where  positive  angles 
are  defined  right  and  up. 

3.8  Scatterer  Motion 

We  will  assume  that  the  yaw,  pitch,  and  roll  motion  of  the  target  is 
confined  to  small  angles.  With  this  assumption  the  target  motion  (excluding 
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translational  motion  of  the  target  c.g.)  will  not  affect  the  location  of  the 
scattering  center  as  far  as  the  ability  of  the  radar  to  resolve  or  measure 
it.  We  do  not  have  to  take  this  motion  into  account  when  computing  the  lo¬ 
cation  of  the  scattering  center. 

However,  the  radar  is  sensitive  to  motion  of  the  scattering  center.  For 
small  angles,  the  motion  of  the  scattering  center  in  the  target  coordinate 
system  is  given  by 


\  ■  *k<l  -  *kR 

(48) 

•j’k’V-V 

(49) 

• 

\  ■  V  -  V 

(50) 

where  P,Q,  and  R  are  the  roll,  pitch,  and  yaw  body  rates  defined  in  Section  3.1. 
The  differential  range  rate  of  the  scatterer  induced  by  target  motion  is  given 
•  by  the  derivative  of  (44)  as 

Ark  “  “Aak  “  +  ykWy(I)  +  YZ(I)]  (51) 

^  This  quantity  is  to  be  added  to  the  range  rate  of  the  target  c.g.,  which  is 

the  derivative  of  (21)  as 


r  -  (<x  -  x  )(x  -x  )  +  (y  -  y  ) (y  -y  )  +  (z  -z  )(z  -z  )]/r 
o  m  o  m  o  m  o  m  o  mono  o 


(52) 


3.9  Effective  Radiated  Power  for  the  RFSS 

The  problem  is  to  simulate  a  target  on  the  RFSS  array  so  that  the  power 
received  by  the  missile  radar  under  test  will  be  equivalent  to  what  would  be 
received  by  the  actual  target.  From  the  radar  range  equation  for  a  reference 
point  scatterer  we  can  write 


PR- 


P tg2x2<j 
(47r)2r4 


(53) 


where  PT  is  the  peak  transmit  power,  G  is  the  one-way  power  gain  of  the  an¬ 
tenna,  X  is  the  wavelength,  r  is  the  range,  and  o  is  the  radar  cross  section 
(RCS)  of  the  point  scatterer.  We  assume  that  the  antenna  boresigbt  is  pointing 
at  the  scatterer. 
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4.  GENERATING  THE  MODULATION  SIGNALS 

In  Section  3  a  ray  has  been  computed  to  each  scatterer  of  the  extended 
target  model  at  one  instant  of  time.  The  information  comput  ch  scat¬ 

terer  is  range,  range  rate,  two  angles,  and  an  amplitude.  In  this  section  we 
will  generate  the  modulation  signals  on  the  basis  of  the  phasor  summation  of 
the  component  rays.  The  procedure  is  based  on  a  constant-Doppler  assumption 
for  each  scatterer  during  the  interval  between  updates  (the  time  between  cal¬ 
culations  of  range,  range  rate,  and  amplitude  for  each  scatterer). 

4.1  Delay  and  Doppler  Coefficients  of  Each  Scatterer 

Regardless  of  the  approach  used  to  generate  the  modulation  signals ,  one 
must  begin  with  the  location  of  each  scatterer  in  radar  coordinates.  In  the 
notation  of  Section  3  the  round-trip  delay  of  the  kc  scatterer  is  given  by 


Tk  -  2(ro  +  Ark)/c 


(58) 


where  r  is  the  range  to  the  target  c.g. ,  Ar,  is  the  differential  range  of 

th°  ^ 

the  k  scatterer,  and  c  is  the  propagation  velocity.  Even  though  the  scat- 
terers  may  be  moving  between  updates  we  will  assume  that  the  radar  is  incapable 
of  measuring  this  motion  on  the  basis  of  its  range  resolution  capability  so 
that  (58)  can  be  assumed  to  be  a  constant  that  will  be  updated  periodically. 

The  Doppler  frequency  of  the  k  scatterer  is  given  similarly  by 


vk  -  -2(ro  +  Ark)A 


(59) 


where  r  is  the  range  rate  of  the  target  c.g.,  Ar  is  the  differential  Doppler 

°th 

of  the  k  scatterer,  and  X  is  the  wavelength.  Even  though  the  scatterers 
may  be  accelerating  between  updates  we  will  also  assume  that  the  radar  is 
Incapable  of  measuring  this  motion  on  the  basis  of  its  Doppler  resolution 
capability.  Thus  (59)  can  be  assumed  to  also  be  a  constant  that  will  be 
update*,  periodically. 
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For  the  moment,  let  us  Ignore  the  discussion  in  Section  2  on  the  medium- 
range  target  model.  Instead,  let  us  assume  that  we  can  apply  an  arbitrary 
delay  and  angular  position  in  space  to  every  scatterer  in  the  target  model 
(such  an  implementation  would  be  straightforward  in  a  nonreal-time  digital 
simulation).  Then  if  A,  represents  the  amplitude  of  the  scatterer,  the 
instantaneous  phasor  assigned  to  the  k  scatterer  is  given  by 

\  eJ2ir(Tk  +  Vkt)  (60) 

2  2 
Mote  that  if  A^  has  the  dimensions  of  RCS,  then  |v^|  does  also. 

An  implementation  of  this  "exact"  approach  requires  a  separate  signal 
generator  for  each  scatterer.  The  delay  and  Doppler  of  each  signal  are  given 
by  (58)  and  (59) ,  and  the  angular  position  in  space  of  the  kc  scatterer  rel¬ 
ative  to  the  target  c.g.  is  given  by  (-Ab^/r^  -  Lc^/tq),  where  Ab^  and  Ac^ 
are  defined  in  Section  3. 

Since  this  approach  is  not  capable  of  being  implemented  on  the  RFSS,  at 
least  when  the  number  of  scatterers  is  large,  we  will  now  concentrate  on  one 
method  that  is  practical. 

4.3  The  Tapped-Delay  Line  Approach 

The  principal  problem  with  the  exact  approach  is  that  the  scatterers 
can  occur  at  arbitrary  delays.  It  is  a  problem  not  only  for  generating  ana¬ 
log  signals,  but  also  in  an  efficient  implementation  of  a  digital  simulation. 
We  solve  this  problem  by  constraining  the  scatterers  to  occur  at  one  of  a 
set  of  discrete  times  in  delay  that  are  uniformly  spaced.  But  in  solving  one 
problem  we  create  another. 

If  we  move  a  scatterer  in  delay  from  its  original  value  in  (58)  to  some 
other  value,  the  radar  may  be  able  to  sense  (measure)  the  change,  especially 
if  the  sample  spacing  of  the  discrete  delays  is  not  small  relative  to  the 
measurement  accuracy  of  the  radar  (as  will  be  the  case  in  any  practical  im¬ 
plementation).  Thus  one  cannot  simply  move  the  scatterer  without  changing 
the  resulting  response  in  the  radar.  However,  one  can  implement  a  resampling 
technique  that  is  described  in  Reference  1  (Section  5). 


This  technique  Is  best  viewed  with  respect  to  Figure  10,  where  four  taps 
of  a  tapped  delay  line  are  located  in  delay  about  the  delay  of  the  scatterer 
of  interest,  two  taps  on  each  side  of  the  scatterer.  We  will  create  new  sig¬ 
nals  at  these  taps  so  that  if  the  range  gate  in  the  receiver  is  centered  at 
the  delay  of  any  one  of  these  taps,  the  receiver  will  measure  the  same  response 
as  it  would  for  the  original  scatterer.  To  find  these  four  signals,  which  we 
designate  as  i*l,  ...,  4,  let  us  define  the  range  gate  response  amplitude, 
X(t),  to  be  a  function  of  delay  normalized  to  the  tap  spacing.  Then  if  p  is 
the  differential  delay  between  the  second  tap  and  the  original  scatterer  we 
can  write  four  equations  as 


X(Hp) 

-  uijX(O)  +  “2XU)  +  + 

(61) 

X(p) 

-  o^XC-l)  +  u2X<°>  +  “3X(D  +  w4X(2) 

(62) 

X(l-P) 

-  u)jX(-2)  +  o>2x(-l)  +  w3x(0)  +  u4x(l) 

(63) 

X<2-p) 

-  di1X(-3)  +  w2X(-2)  +  w3x(-l)  +  <*>4x(°) 

(64) 

In  other  words,  the  left  side  of  (61)  represents  the  amplitude  weighting  that 
would  be  placed  on  the  original  scatterer  if  the  range  gate  were  centered  on 
the  first  tap,  while  the  right  side  represents  the  composite  weighting  that 
would  be  applied  to  the  four  signals  originating  at  the  taps.  We  can  solve 


{u^,  ...,  u)^}  «  tap  weights 


r 

1^2  f  w 


original  scatterer 


A 

At-1  ["*- 


Delay 


Figure  10.  Creation  of  Signals  at  Four  Taps  of  Tapped-Delay  Line 


this  set  of  equations  for  through  and  the  solution  will  apply  to  a 
specific  value  of  p.  The  solution  obtained  in  (61)  through  (64)  is  valid  as 
long  as  the  tap  spacing  is  about  5QZ  of  the  range  gate  spacing  in  the  radar. 
We  also  assume  that  the  range  gate  response  drops  off  sufficiently  fast  so 
that  no  more  than  four  taps  (and  four  equations)  need  be  evaluated. 

4.4  Nation  of  Signals  at  Each  Tan 

At  each  tap  three  signals  are  formed ,  where  each  signal  is  a  weighted 
summation  of  the  phasors  defined  by  (60) .  The  first  signal  for  the  nC^ 
tap  is  the  scintillation  signal  given  by 


v  r  _  ,  v  -j4TT(i-n+l+p)AR_A 

»<«)  -  I  'kVw<f)*  ^ 


(65) 


where  the  summation  is  over  the  set  of  scatterers ,  AR^  is  the  tap  spacing  in 


range  units  (ARj.  -  cAtt/2  where  Att  is  the  delay  between  taps) ,  l  is  the  tap 
number  for  the  first  of  four  that  surround  the  kth  scatterer,  and  p  is  the 


fractional  distance  the  1c  scatterer  is  from  the  second  of  four  taps  that 
surround  the  scatterer.  The  tap  weights  are  nonzero  only  for  i-l,...,4. 
The  phase  factor  in  (65)  is  to  compensate  for  the  fact  that  the  original 


scatterer  is  replaced  by  other  components  at  different  ranges. 

The  second  signal  for  the  nth  tap  is  the  angle-weighted  response  in  the 


horizontal  direction  (oriented  to  the  RFSS  chamber,  or  the  B-direction  in  the 
notation  of  Section  3.4).  It  is  given  by 


v°>  -  f  i  *vkW’>*-wl'“+l+p>iVX 


(66) 


■•o  k 


The  third  signal  is  the  angle-weighted  response  in  the  vertical  direction 
(oriented  to  the  RFSS  chamber)  and  it  is  given  by 
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Vn> 


f  s* 


,ckVk“n-i,+l 


(p  )  e-i4ir  ( Jl-n+i+p  )  A VX 


(67) 


4.5  Computation  of  Glint  Offsets  at  Each  Tap 

The  angular  glint  offsets  for  each  tap  are  computed  as 


A^n)  -  -Re{Vg  (n)  /V  (n)  } 

(68) 

A^n)  -  -Re{Vc(n)/V(n)} 

(69) 

where  we  define  the  azimuth  glint  offset 
hand  horizontal  direction  (referenced  to 
B-direction) ,  and  elevation  glint  offset 
negative  C -direction) . 


A^fa)  to  be  positive  in  the  right- 
the  RFSS  chamber,  or  in  the  negative 
A_  (n)  to  be  positive  up  (or  in  the 

LL 


f 
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5.  SIMULATION  ARCHITECTURE 


The  computer  resources  that  are  available  at  the  KFSS  to  generate  ex¬ 
tended  target  signals  consist  of  a  Datacraft/1  minicomputer  that  acts  as  a 
host  to  the  AP120B  array  processor  built  by  Floating  Point  Systems.  The 
input  to  the  host  computer  consists  of  the  state  of  the  engagement  geometry 
and  the  output  of  the  AP120B  consists  of  the  Doppler  modulation  signal  and 
glint  offsets  at  each  tap  (V(n),  A^ (n) ,  A^Cn) ) .  In  deciding  which  opera¬ 
tions  are  to  be  assigned  to  the  two  processors,  the  following  factors  were 
applied: 

1)  the  array  processor  works  best  on-  arrays  where  the  number  of 
operations  on  each  sample  is  few; 

2)  the  program  construction  is  costly  for  the  array  processor 
so  that  program  should  be  fairly  stable; 

3)  modifications  are  best  made  in  the  host  computer;  and 

4)  there  should  be  relatively  little  traffic  over  the  real-time 
interface  between  the  host  computer  and  the  array  processor. 

With  these  -factors  in  mind,  the  assignment  of  the  various  processing  steps 
to  generate  extended  target  signals  reduces  naturally  to: 

The  host  computer 

1.  transform  engagement  geometry 

* 

2 .  compute  visibility  and  amplitude  of  scatterers 
The  array  processor 

1.  compute  which  taps  are  affected  by  each  scatterer 

2.  compute  tap  weights 

3.  compute  and  sum  phasors  at  output  sample  rate 

4.  compute  the  Doppler  modulation  signal  and  glint  offsets  at  each 

tap 


In  some  cases  the  computation  for  the  amplitude  of  the  scatterers  might 
be  too  long  for  implementation  in  the  host  computer.  In  such  a  case,  this 
f  computation  can  be  transferred  to  the  AF120B. 


> 


The  real-time  interface  between  the  host  and  the  AP120B  consists  of 
tQ  and  r^»  the  range  to  the  target  c.g.  and  the  range  to  the  first  tap; 

Aa^,  AA^»  Ab^>  and  Ac^»  the  ABC  vector  components  for  each  visible  scatterer 
and  the  amplitude  of  the  iC^  scatterer. 

In  Appendix  B  we  describe  and  give  listings  for  an  extended  target  sim¬ 
ulation  program  that  conforms  to  the  above  architecture. 
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APPENDIX  A 


LIMITING  RANGE  FOR  MEDIUM-RANGE 
TARGET  MDDEL 


Introduction 

A  medium-range  target  model  was  developed  In  Reference  1,  consisting 
of  N  point  scatterers,  where  each  scatterer  can  have  an  RCS  that  Is  aspect 
dependent.  The  medium-range  constraint  assumes  that  all  scatterers  on  the 
target  are  In  the  linear  region  of  the  monopulse  receive  beam,  and  all 
scatterers  are  Illuminated  with  a  constant  gain  by  the  transmit  beam.  The 
purpose  of  this  constraint  Is  to  remove  the  sensor  pointing  angles  from  the 
real-time  computation.  In  other  words,  the  signal  that  Is  generated  on  the 
RFSS  array  is  Independent  of  the  sensor  pointing  angles.  At  shorter  ranges 
where  the  transmit  beam  is  no  longer  uniform  across  the  target,  or  where  the 
monopulse  difference  beam  is  not  linear,  the  pointing  angles  of  the  sensor 
beam  must  be  known  so  that  the  variable  weighting  can  be  Implemented  in  the 
real-time  simulation;  moreover,  the  signals  that  would  be  received  on  each 
monopulse  channel  must  be  separately  simulated  and  radiated  into  specific 
points  on  the  receive  beam  so  that  each  channel  receives  the  proper  signal 
and  rejects  the  others. ^ 

The  purpose  of  this  memo  is  to  determine  the  minimum  range  for  the 
applicability  of  the  medium-range  model.  A  simple  Monte-Carlo  simulation 
will  be  used  to  accomplish  this. 

The  Target  Model 

A  statistical  type  target  model  is  assumed.  Two  scatterers  are  separated 
by  an  angle  8^,  in  between  N-2  scatterers  are  placed  at  random.  Thus  the 


^"Design  Requirements  for  Simulating  Realistic  RF  Environment  Signals 
on  the  RFSS,"  MRI  Report  132-44,  by  R.  L.  Mitchell  and  I.  P.  Bottlik,  dated 
23  September  1977. 
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target  consists  of  N  scatterers  that  cover  an  angular  width  of  6^.  All 
scatterers  are  assumed  to  be  of  equal  RCS  on  the  average,  and  each  is 
fluctuated  with  a  Rayleigh  amplitude  and  random  phase. 

The  Antenna  Patterns 

The  stim  channel  two-way  voltage  antenna  pattern  is  assumed  to  be 


Gj.0)  -  1  -  60' 


and  the  two-way  difference  pattern 


Ga(6)  -  k(0-a0  ;  (2) 

Therefore,  if  the  complex  voltage  is  on  the  ith  scatterer  at  an  angle 
0^,  then  the  received  voltages  on  the  two  channels  are 


\  -]C  <i-69i,vi 


\  -  kswi'a9i)vl 


For  the  purpose  of  this  investigation  we  have  assumed 


a  -  1 .70/ 0^ 


6  -  1-37/e^  (6) 

where  0^^  is  the  one-way  half-power  beamwidth.  These  values  are  typical 
of  many  tracking  radars.  The  constant  k  will  factor  out  of  the  problem 
f  later. 


We  assume  Chat  the  boresite  of  the  antenna  Is  pointing  exactly  at  the 
center  of  the  target  (midway  between  the  end  points).  The  estimate  of 
angle  is  assumed  to  be 


9 ACT  "  k  Re{'VV2J 


where  the  subscript  ACT  denotes  the  actual  (assumed)  target ,  in  contrast 
to  an  approximate  one  based  on  the  medium  range  model. 

The  Medium-Range  Model 

A  glint  centroid  will  be  calculated  for  the  target  that  is  based  on 
Gf(9)  *  1  and  G^(9)  ■  8*  Thus  the  composite  signal 


-•  -z* 


will  be  radiated  from  the  angle 


iT,  v 

vZ— 1 


9  -  Re 


Row  if  we  use  the  antenna  patterns  in  (1)  and  (2) ,  and  the  formula  for  the 
estimate  of  the  angle,  we  have 

®4PP  ■  8’  ^4  (10> 

1-0(6  r 

*  e 

Thus  we  will  compare  9^pp  with  9^^  to  determine  where  the  medium  range 
model  breaks  down. 


Results 

In  Tables  A-l  through  A-6  we  show  the  results  of  20  statistical  replica¬ 
tions  of  a  target  consisting  of  N*5  scatterers,  where  the  target  width  varies  from 
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6T/83dB 


.25  to  1.50  (the  glint  angles  0^^  and  0^^  are  designated  as 


ACTUAL  and  APPROX,  each  being  normalized  to  the  half -power  width).  For  a 
target  width  of  252  of  the  beamwidth  (Table  A-l)  the  peak  error  is  .003 
(or  .32  of  the  beamwidth),  which  is  negligible.  For  a  target  width  of 
502  of  the  beamwdith  (Table  A-2)  the  peak  error  is  over  100%  of  the  beamwidth 
(REP  20);  however,  the  actual  glint  for  this  case  is  also  large,  amounting 
to  662  of  the  beamwidth.  In  practice,  we  can  tolerate  a  large  error  if 
the  glint  angle  is  also  large.  It  is  more  important  to  keep  the  errors 
small  when  the  glint  angles  are  small.  Thus  REP  12  in  Table  A-2  represents 
probably  the  most  severe  error,  which  is  2.8%  of  the  beamwidth  when  the 
actual  glint  is  142  of  the  beamwidth. 

If  we  rule  out  those  replications  where  the  actual  glint  is  larger  than 
half  of  the  target  width,  we  can  construct  the  following  table 


Target  Width 

Peak  Error 

.25 

.002 

.50 

.028 

.75 

.076 

1.00 

.184 

are  negligible. 

However,  when  we 

width  1.250-,-  (Table  A-5)  we  observe  several  large  errors,  even  when  the 

JOB 

actual  glint  is  small.  For  example,  on  REP  16  the  actual  glint  is  only 
7.72  of  the  beamwidth,  but  the  error  is  over  2  beamwidths.  Clearly,  the 
model  breaks  down  for  0^,  «  1. 258^2 .  The  actual  point  at  which  the  model 
breaks  down  lies  somewhere  between  8^,  ■  l.OOe^g  and  0T  *  1.250^^. 
Variations  in  the  antenna  patterns  and  formulas  for  measuring  angle  will 
Impact  on  a  precise  determination  of  where  the  model  breaks  down,  but  we 
can  state  conservatively  that  the  model  is  valid  as  long  as  0^  <_ 

In  order  to  test  the  effect  of  the  number  of  scatterers  in  the  model, 
we  repeated  the  previous  simulation  for  N»10.  The  results  are  shown  in 
Tables  A-7  through  A-12.  No  major  discrepancies  are  noted  from  the  previous 
1  conclusions. 
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Table  A-l.  Comparison  of  Actual  Target  with  Medium-Range  Model  (APPROX) 


£ot  9T/93dB 


0.25  (N-5,  all  angles  normalized  to  dOJ_) 


REP 

ACTUAL 

APPROX 

01 FF 

.  SCATTERER  LOG 

ATI  ON.  .  . 

1 

.  051 

.  051 

000 

125 

.  020 

.  072 

.  113 

.  125 

2 

057 

058 

001 

125 

123 

051 

-  012 

.  125 

3 

.  043 

.  043 

.  000 

125 

056 

049 

.  047 

.  125 

4 

074 

073 

.  001 

125 

092 

029 

.  083 

.  125 

5 

175 

178 

-  003 

-  125 

100 

056 

.  021 

.  125 

6 

.  167 

.  165 

001 

125 

104 

.  030 

.  123 

.  125 

7 

.  168 

.  169 

.  001 

125 

.  048 

.  109 

.  120 

.  125 

8 

.  002 

.  002 

000 

125 

092 

072 

.  091 

.  125 

9 

.  015 

.  015 

.  000 

125 

-  120 

.  010 

.  080 

.  125 

10 

-  018 

-  019 

000 

125 

— .  046 

.  C66 

.  110 

.  125 

11 

.  041 

.  041 

.  000 

125 

-  053 

.  021 

.  101 

.  125 

12 

.  095 

.  096 

.  001 

125 

026 

001 

.  097 

.  125 

13 

-  045 

047 

002 

125 

.  024 

.  093 

.  108 

.  125 

14 

.  035 

.  035 

.  000 

125 

-.  073 

009 

.  097 

.  125 

15 

069 

071 

002 

125 

013 

Oil 

.  032 

.  125 

16 

035 

034 

.  001 

125 

066 

.  023 

.  066 

.  125 

17 

.  077 

.  077 

000 

125 

014 

001 

.  105 

.  125 

18 

109 

110 

001 

125 

087 

072 

.  083 

.  125 

19 

.  113 

.  114 

.  000 

-  125 

.  044 

.  087 

.  107 

.  125 

20 

006 

006 

.  000 

125 

062 

055 

035 

.  125 

Table  A-2.  Comparison  of  Actual  Target  with  Medium-Range  Model  (APPROX) 


£ot  9T/93dB 


0.50  (N**5,  all  angles  normalized  to 


REP 

ACTUAL 

APPROX 

D1FF 

.  SCATTERER  LOCATION.  .  . 

1 

096 

-  102 

006 

250 

203 

094 

061 

.  250 

2 

042 

054 

Oil 

-  250 

.  005 

.  075 

.  083 

.  250 

3 

220 

208 

.  012 

250 

192 

-  156 

057 

.  250 

4 

.  122 

.  126 

.  004 

250 

.  105 

.  Ill 

.  176 

.  250 

5 

940 

-1.  774 

-  835 

250 

242 

214 

.  071 

.  250 

6 

459 

445 

.  014 

250 

056 

039 

.  227 

.  250 

7 

172 

171 

.  001 

250 

212 

.  052 

.  199 

.  250 

8 

.  116 

.  110 

006 

250 

.  036 

.  038 

.  204 

.  250 

9 

263 

277 

014 

250 

189 

.  045 

.  142 

.  250 

10 

.  057 

.  063 

.  006 

250 

133 

Ill 

024 

.  250 

11 

.  047 

.  044 

-.  003 

250 

181 

.  062 

.  152 

.  250 

12 

.  141 

.  113 

-.  028 

-.  250 

.  C08 

.  067 

.  120 

.  250 

/  13 

.  088 

.  087 

001 

250 

.  C04 

.  101 

.  124 

.  250 

14 

.  031 

.  036 

.  004 

250 

096 

007 

.  031 

.  250 

15 

— .  666 

-.  400 

.  266 

250 

229 

.  150 

.  188 

.  250 

16 

.  381 

.  364 

017 

250 

.  177 

.  177 

.  199 

.  250 

17 

305 

234 

.  021 

250 

149 

.  001 

.  094 

.  250 

18 

.  032 

.  040 

.  008 

250 

158 

101 

.  179 

.  250 

19 

.  075 

.  076 

.  001 

250 

189 

141 

.  156 

.  250 

SO 

.  662 

1.  777 

1.  115 

-  250 

.  018 

.  039 

.  227 

.  250 
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Table  A-3.  Comparison  of  Actual  Target  with  Medium-Range  Model  (APPROX) 


for  eT/63dB  -  1.25  (N-5, 


Rt^ 

ACTUAL 

AFrROX 

DIFF 

1 

315 

332 

017 

— 

2 

217 

079 

.  133 

- 

3 

1.  017 

-2.  253 

-3.  271 

— 

4 

.  229 

.  245 

.  016 

— 

5 

.  038 

-.  007 

045 

— 

6 

.  459 

.  257 

-.  202 

- 

7 

118 

032 

.  036 

— 

a 

394 

-1.  867 

-1.  473 

— 

9 

-.  054 

039 

036 

— 

10 

.  535 

3.  309 

2.  774 

— 

11 

.  496 

.  241 

255 

— 

12 

.  222 

131 

352 

— 

13 

429 

1.  867 

2.  296 

— 

14 

051 

.  216 

.  267 

• 

15 

-  005 

•  152 

.  147 

— 

16 

.  077 

2.  433 

2.  406 

— 

17 

-1.  712 

279 

1.  433 

— 

18 

400 

453 

053 

— 

19 

.  010 

.  018 

.  007 

— 

20 

260 

334 

124 

—  # 

all  angles  normalized  to  0OJ_) 

oaB 

. SCATTERER  LOCATION.  . 


625 

549 

540 

369 

.  625 

625 

453 

-.  304 

.  399 

.  625 

625 

553 

-  303 

230 

.  625 

625 

549 

.  170 

.  410 

.  625 

625 

240 

.  028 

.  241 

.  625 

625 

227 

.  348 

.  516 

.  625 

625 

220 

062 

.  no 

.  625 

625 

573 

-.  105 

.  085 

.  625 

625 

524 

.  033 

.  192 

.  625 

625 

.  136 

.  188 

.  561 

.  625 

625 

.  158 

.  508 

.  553 

.  625 

625 

211 

.  286 

.  466 

.  625 

625 

363 

135 

.  105 

.  625 

625 

211 

098 

.  *39 

.  625 

625 

253 

173 

.  082 

.  625 

625 

205 

.  024 

.  550 

.  625 

625 

150 

.  068 

.  156 

.  625 

625 

377 

302 

271 

.  625 

625 

032 

.  028 

.  605 

.  625 

625 

534 

216 

073 

.  625 

Table  A-4.  Comparison  of  Actual  Target  with  Medium-Range  Model  (APPROX) 
for  *  1.50  (N»5,  all  angles  normalized  to  03dfi) 


REP 

ACTUAL 

APFROX 

DIFF 

1 

.  476 

155 

-.  631 

2 

.  139 

236 

-.  375 

3 

091 

14.  613 

14.  704 

4 

093 

2.  506 

2.  598 

5 

138 

-.  129 

.  009 

6 

.  034 

-.  245 

-.  280 

7 

064 

.  343 

.  407 

8 

-1.  985 

-.  294 

1.  691 

9 

.  090 

026 

-.  117 

10 

.  229 

016 

-.  245 

11 

351 

1.  839 

2.  190 

12 

160 

.  238 

.  443 

13 

-.  396 

-2.  593 

-2.  198 

14 

.  261 

1.  732 

1.  521 

15 

.  155 

1.  731 

1.  626 

16 

.  982 

-8.  846 

-9.  829 

17 

028 

1.  826 

1.  354 

18 

.  023 

221 

244 

19 

-.  005 

.  433 

.  438 

20 

-.  352 

301 

.  051 

SCATTERER  LOCATION 


750 

-.  294 

.  204 

.  466 

.  750 

-.  750 

634 

299 

.  029 

.  750 

750 

-.  542 

-.  200 

.  463 

.  750 

750 

518 

-.  348 

.  004 

.  750 

750 

544 

352 

030 

.  750 

-.  750 

-.  637 

-  264 

.  575 

.  750 

-.  750 

296 

.  097 

.  380 

.  750 

750 

-.  510 

.  558 

.  618 

.  750 

-.  750 

264 

-.  009 

.  398 

.  750 

-.  750 

-.  309 

.  302 

.  347 

.  750 

-.  750 

-.  471 

-.  379 

.  746 

.  750 

-.  750 

302 

.  676 

.  683 

.  750 

750 

-.  548 

.  535 

.  598 

.  750 

-.  750 

641 

-.  018 

.  679 

.  750 

-.  750 

686 

.  555 

.  704 

.  750 

-.  750 

-.  691 

147 

.  341 

.  750 

750 

— .  266 

232 

.  749 

.  750 

-.  750 

-.  743 

295 

.  678 

.  750 

-.  750 

-.  703 

016 

.  748 

.  750 

750 

479 

.  386 

.  740 

.  750 
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Table  A-5.  Comparison  of  Actual  Target  with  Medium-Range  Model  (APPROX) 
for  ®^®3dB  *  (N»5,  all  angles  normalized  to  e^g) 


I 


REP  ACTUAL  APPROX  DIFF 


SCATTERER  LOCATION 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


981 

-1.  775 

794 

375 

470 

454 

.  016 

375 

.  186 

.  171 

015 

375 

.  258 

.  284 

.  025 

375 

.  317 

.  331 

.  013 

375 

.  345 

.  346 

.  001 

375 

.  675 

.  443 

232 

375 

.  154 

.  181 

.  027 

-  375 

.  223 

.  253 

.  030 

375 

168 

iso 

013 

375 

.  235 

.  311 

.  076 

375 

.  201 

.  214 

.  013 

375 

.  163 

.  172 

.  010 

375 

-  023 

013 

.  009 

375 

103 

113 

010 

375 

.  435 

.  438 

.  003 

375 

.  040 

.  045 

.  004 

375 

1.  826 

8.  418 

6.  592 

375 

.  114 

.  151 

.  036 

375 

.  261 

.  291 

.  030 

375 

103 

.  108 

.  173 

.  375 

037 

.  004 

.  290 

.  375“ 

278 

.  198 

.  286 

.  375 

319 

.  199 

.  312 

.  375 

.  005 

.  016 

.  199 

.  375 

.  314 

.  339 

.  343 

.  375 

134 

.  188 

.  228 

.  375 

112 

.  252 

.  309 

.  375 

031 

.  032 

.  226 

.  375 

154 

118 

.  276 

.  375 

339 

233 

.  150 

.  375 

325 

311 

.  180 

.  375 

295 

006 

.  022 

.  375 

348 

236 

096 

.  375 

.  092 

.  146 

.  315 

.  375 

.  072 

.  334 

.  342 

.  375 

197 

122 

.  018 

.  375 

064 

.  067 

.  096 

.  375 

232 

041 

.  128 

.  375 

.  017 

.  077 

.  079 

.  375 

Table  A-6. 


Comparison  of  Actual  Target  with  Medium-Range  Model  (APPROX) 
for  "  1.00  (N**5,  all  angles  normalized  to  0^g) 


REP 

ACTUAL 

AFPRDX 

DIFF 

. SCATTERER  LOCATION.  .  . 

1 

170 

128 

.  042 

500 

352 

230 

.  350 

.  500 

2 

.  129 

.  174 

.  044 

500 

194 

122 

.  005 

.  500 

3 

.  283 

.  336 

•  043 

500 

319 

237 

.  316 

-  500 

4 

.  172 

.  356 

.  ie4 

500 

497 

198 

.  481 

.  500 

5 

289 

277 

•  012 

500 

448 

282 

207 

.  500 

6 

247 

300 

053 

500 

395 

.  024 

.  037 

.  500 

7 

.  230 

.  182 

048 

500 

454 

.  090 

.  247 

.  500 

8 

438 

437 

.  001 

500 

434 

424 

091 

.  500 

9 

1 14 

252 

138 

500 

096 

058 

.  295 

.  500 

10 

1.  003 

1.  775 

.  772 

500 

088 

061 

.  259 

.  500 

1 1 

.  221 

.  216 

005 

500 

330 

.  057 

.  112 

.  500 

12 

.  569 

.  071 

499 

500 

299 

.  154 

.  227 

.  500 

13 

018 

004 

.  015 

500 

410 

150 

.  438 

.  500 

14 

.  352 

.  221 

130 

500 

221 

052 

.  492 

.  500 

i  9 

.  149 

.  098 

050 

500 

346 

-  313 

.  180 

.  500 

16 

397 

326 

.  071 

500 

377 

.  295 

.  455 

.  500 

17 

341 

369 

028 

500 

476 

078 

.  178 

.  500 

18 

— .  1 16 

097 

.  020 

500 

466 

226 

.  081 

.  500 

19 

.  264 

.  363 

.  099 

500 

336 

294 

.  098 

.  500 

20 

.  398 

.  416 

.  018 

500 

161 

152 

.  309 

.  500 

A 
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Table  A-7.  Comparison  of  Actual  Target  vith  Medium-Range  Model  (APPROX) 
for  ®j/®3<jb  ■  -25  (N«10,  all  angles  normalized  to  Q^^) 


REP 

ACTUAL 

APPROX 

DIFF 

1 

.  014 

.  014 

.  001 

2 

.  033 

.  033 

.  000 

3 

029 

029 

.  000 

4 

.  042 

.  042 

000 

5 

.  032 

.  031 

001 

6 

007 

007 

.  000 

7 

000 

000 

000 

8 

014 

C06 

.  007 

9 

.  010 

.  010 

000 

10 

092 

096 

003 

11 

023 

023 

.  000 

12 

.  139 

.  138 

001 

13 

.  046 

.  047 

.  001 

14 

.  042 

.  042 

.  000 

13 

055 

057 

001 

16 

.  018 

.  019 

.  001 

17 

.  003 

.  003 

000 

18 

118 

117 

.  000 

19 

136 

139 

-  003 

20 

056 

058 

002 

Table  A-8.  Comparison  of  Actual  Target  with  Medium-Range  Model  (APPROX) 
for  ®T/®3 £3  *  *50  (N®10,  all  angles  normalized  to  e^g) 


REP 

ACTUAL 

APPROX 

DIFF 

1 

.  049 

.  060 

.  on 

2 

.  040 

.  048 

.  008 

3 

.  167 

.  17B 

.  Oil 

4 

.  242 

.  272 

.  030 

5 

.  177 

.  191 

.  014 

6 

.  042 

.  042 

000 

7 

110 

108 

.  002 

8 

038 

045 

006 

9 

.  032 

.  037 

.  005 

10 

1.  030 

1.  792 

.  762 

11 

039 

049 

009 

12 

.  Ill 

.  122 

.  010 

13 

136 

147 

.  009 

14 

.  161 

.  167 

.  006 

15 

052 

055 

002 

16 

206 

218 

012 

17 

.  028 

.  027 

001 

18 

217 

208 

.  009 

19 

095 

098 

004 

20 

121 

103 

.  018 
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Table  A-9.  Comparison  of  Actual  Target  with  Medium-Range  Model  (APPROX) 
for  ■  >75  (N-10,  all  angles  normalized  to  Q^g) 


REP 

ACTUAL 

APPROX 

DIFF 

1 

049 

043 

.  006 

2 

.  017 

.  000 

016 

3 

.  022 

.  026 

.  004 

4 

141 

142 

001 

5 

.  542 

.  451 

092 

6 

083 

-.111 

028 

/ 

.  021 

.  023 

.  002 

8 

.  226 

.  108 

118 

9 

1.  156 

1.  825 

.  669 

10 

.  188 

.  224 

.  036 

11 

043 

043 

.  000 

12 

.  282 

.  356 

.  074 

13 

.  055 

.  053 

001 

14 

.  119 

.  093 

026 

IS 

.  192 

.  228 

.  036 

16 

232 

303 

072 

17 

.  587 

.  036 

551 

18 

.  247 

.  291 

.  043 

19 

.  036 

.  031 

005 

20 

309 

299 

.  010 

Table  A-10. 


Comparison  of  Actual  Target  with  Medium-Range  Model  (APPROX) 


for  ®T/®3dB 


1.00  (N»10,  all  angles  normalized  to  e^g) 


ACTUAL 

APPROX 

DIFF 

045 

096 

051 

.  048 

.  068 

.  020 

013 

.  003 

•  016 

.  681 

2.  162 

1.  481 

426 

442 

016 

.  327 

.  364 

.  037 

045 

-  128 

083 

.  129 

.  232 

.  103 

411 

458 

047 

.  055 

.  077 

.  022 

173 

073 

.  100 

.  143 

.  124 

019 

.  292 

.  391 

.  099 

109 

154 

045 

.  224 

.  314 

.  090 

201 

227 

025 

.  462 

.  459 

003 

346 

366 

020 

.  203 

.  414 

•  211 

036 

021 

.  015 
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Table  A-ll.  Comparison  of  Actual  Target  with  Medium-Range  Model  (APPROX) 
for  *  1*25  (N-10,  all  angles  normalized  to  9,JB) 


3dB 


REP  ACTUAL  APPROX  DIFF 


Table  A- 12. 


1 

220 

220  001 

2 

.  368 

■  402  .  034 

3 

-  271 

365  094 

4 

.  077 

•  133  .  061 

5 

.  187 

258  445 

6 

184 

-■  455  271 

7 

.  027 

1.922  1.895 

8 

006 

-.  082  076 

9 

284 

356  072 

10 

155 

233  -  078 

11 

-.  811 

-2.  001  -1.  190 

12 

229 

277  -.  048 

13 

.  248 

•  216  -.032 

14 

.  041 

025  065 

IS 

363 

-■  431  069 

16 

070 

-•  177  107 

17 

675 

-3.  174  -2.  499 

18 

136 

•  441  .  577 

19 

058 

167  108 

20 

.  375 

•  357  018 

son  of  Actual  Target  with  Medium-Range  Model 

93dB  “ 

1.50  (N“1G,  all  angles  normalized  to  1 

REP 

ACTUAL 

APPROX  DIFF 

1 

.  153 

•  418  .  265 

2 

.  053 

•  156  .  103 

3 

.  800 

•  431  -.369 

4 

.  242 

•  456  .  214 

5 

165 

■  170  .  335 

6 

.  068 

269  337 

7 

222 

“•  439  -.217 

8 

.  155 

•  452  .  297 

9 

.  048 

“■  254  -.  302 

10 

.  017 

•  284  .  268 

11 

.  863 

■  095  -.  768 

12 

.  006 

-2.  603  -2.  609 

13 

.  179 

-■  085  -.  263 

14 

.  006 

“■  173  -.  179 

15 

003 

•  145  .  148 

16 

305 

-1.  774  -1.  470 

17 

.  019 

“■  447  -.  466 

18 

.  057 

•  213  .  156 

19 

044 

-•  330  -.  336 

20 

.  140 

. 204  .  064 

3dB'> 


APPENDIX  B 


FORTRAN  PROGRAM  FOR  GENERATING  REAL-TIME 
EXTENDED  TARGET  SIGNALS 

The  Fortran  program  described  here  generates  the  Doppler  modulation 
signal  and  glint  offsets  at  each  tap  of  a  tapped-delay  line  for  an  extended 
target  composed  of  a  set  of  discrete  scatterers.  It  represents  the  latest 
version  delivered  to  MIRADCOM  on  20  June  1978,  and  only  minor  corrections 
in  comment  statements  are  made  in  the  following  listings,  with  the  exception 
of  subroutine  XFORM  (where  the  sign  convention  on  all  angles  was  reversed). 

There  are  two  principal  subroutines  in  this  program,  TARGEO  which  com¬ 
putes  the  amplitude  and  geometrical  information  for  each  scatterer,  and 
TARGDM  which  computes  the  Doppler  modulation  and  glint  offsets  for  each  tap. 
This  architecture  assumes  that  TARGEO  will  be  installed  in  the  host  computer 
and  TARGDM  in  the  AP120B  array  processor.  Block  diagrams  for  these  two  sub¬ 
routines  are  sketched  in  Figures  B-l  and  B-2.  The  assignment  of  the  other 
subroutines  is  as  follows: 

Host  Computer 

MAIN  -  main  or  driver  program 

ETGE0  -  updates  and  transforms  engagement  geometry 

XFORM  -  transforms  inertial  coordinates  to  target  coordinates 

SCTAMP  -  computes  amplitude  of  scatterer 

AP120B 

ETGDM  -  generates  Doppler  modulation  and  glint  offsets 
TAPUTS  -  compute  tap  weights  (table  lookup) 

Initialization  (Host) 

TAPSET  -  generates  tap  weight  table  (for  TAPWTS) 

CHI  -  range  gate  response  (used  by  TAPSET) 

SIMQ  -  simultaneous  equation  solution  (used  by  TAPSET) 

DATAIN  -  reads  target  scattering  data  from  cards 
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*IF  VISIBLE  N„  »  NUMBER  OF 

S  VISIBLE 
SCAITERERS 


Figure  B-l.  BLOCK  DIAGRAM  FOR  SUBROUTINE  TARGEO 

(coopuCes  amplitude  &  geometrical  info  for  each 
scatterer) 


Utility 

XMIT  -  moves  data  ^ 

S INC OS  -  computes  sine/cos 

A  sample  test  program  Is  Included  as  MAIN  in  the  following  listings t  and  an 

alternate  Doppler  modulation  and  glint  offset  subroutine  (ETGD1)  is  Included 

that  Is  based  on  a  single  tap  (no  range  extent) .  4 


NOTE:  The  programs  listed  here  for  eventual  installation  on  the  AP120B  are 
written  in  FORTRAN;  they  must  be  converted  to  the  appropriate  language  on  the 
SP120B  in  order  to  run  in  real  time. 


a  a  a  a  a  a  a  a  a  a  a  o  o  n  o  a  a  a  o  a  a  o  o  a 
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PROGRAM  MAIN  < INPUT, OUTPUT,  TAPE5=INPUT. TAPE6=0UTPUT> 

THIS  IS  A  SAMPLE  MAIN  PROGRAM  FOR  TEST  PURPOSES  ONLY. 

THIS  EXTENDED-TARGET  SIMULATION  PACKAGE  HAS  BEEN  PREPARED  BY 
RL  MITCHELL  OF  MARK  RESOURCES,  INC  (213-S22-4955),  UNDER  CONTRACT  TO 
MIRADCOM.  IT  IS  WRITTEN  WITH  THE  INTENTION  THAT  IT  WILL  BE  MADE  PART 
OF  A  REAL-TIME  SIMULATION  PROGRAM,  ALTHOUGH  SOME  OF  THE  CODE  IS  NOT 
WRITTEN  IN  COMPLETELY  OPTIMUM  FORM  (IT  IS  MORE  EASILY  UNDERSTOOD  THIS 
WAY,  AND  THE  REVISIONS  ARE  EASILY  MADE). 

ALL  ARRAYS  IN  COMMON  SHOULD  BE  DIMENSIONED  IN  THE  MAIN  PROGRAM. 

SEE  THE  SUBROUTINES  FOR  A  DEFINITION  OF  THE  VARIABLES. 

RULES  FOR  DIMENSIONING  ARRAYS . 

X,  Y,  Z . .  NSC  AT 

AMP, DA. DB, DC, DAD.  .  .  .  NSCAT  (MAYBE  SMALLER) 

TWARAY . 4*NARAY 

VR,  VI.  DAZ.  DEL . NTAP 


COMMON  /Tl/  XO,  YO,  ZO.  XOD,  YOD,  ZOD,  PSI,  THETA,  PHI.  BP.  BQ,  BR 

COMMON  /T2/  CPSI. SPSI,  CTHETA,  STHETA,  CPHI, SPHI 

COMMON  /T3/  XM.  YM,  ZM,  XMD,  YMD,  ZMD 

COMMON  /T4/  NSCAT. ST.  AMPMIN,  X(20),  Y(20), Z(20) 

COMMON  /T5/  NTAP, DRTAP,  DRGATE,  XL,  PTGDSQ 
COMMON  /T6/  NAR AY, TWARAY (404) 

COMMON  /Tl  /  NS,  RO,  R 1 ,  ROD, AMP ( 20) ,  DA ( 20) , DB ( 20 ) , DC ( 20) , DAD ( 20) 
COMMON  /T8/  PEFF, VR (8) , VI (8) ,  DAZ ( 8) ,  DEL (8) 


DEFINE  VARIABLES. 


DATA  LR.  LW/5.  6/ 

DATA  NTAP/8/,  DRTAP/30.  /,  DRCATE/60.  /,  XL/.  02/,  AMPMIN/ 1.  E-10/. 
I  PTGDSQ/ 1.  / 

DATA  DTI ME/ 1.  / 

DATA  NARAY/I01/ 

DATA  ST/O.  / 


READ 

READ 

READ 

READ 

READ 

READ 

READ 


(LR, 

(LR, 

(LR, 

(LR, 

(LR, 

(LR, 

(LR, 


XO,  YO,  ZO 
XOD, YOD, ZOD 
100)  PSI,  THETA,  PHI 
100)  BP,  BQ,  3R 
XM.  YM,  ZM 
XMD,  YMD,  ZMD 


100) 

100) 


100) 

100) 


101)  NSCAT 


J 


ooo 
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< 


READ 


<LR, 100) 


4 


(X(K),  Y(K).  Z(K).  K=l, NSCAT) 


WRITE  (LW,  200) 
WRITE  (LW,  201 ) 
WRITE  (LW, 202) 
WRITE  (LW, 203) 
WRITE  (LW, 204) 
WRITE  (LW,  205) 
WRITE  (LW, 206) 
WRITE  (LW, 207) 


XO,  YO,  20 
XOD,  YOD,  20D 
PSI, THETA, PHI 

3P,  BQ, BR  '  4 

XM,  YM,  ZM 
XMD. YMD, ZMD 
NSCAT  • 

(X(K), Y(K), Z(K).  K-l.  NSCAT) 


SUBROUTINES  TO  BE  CALLED  FROM  MAIN  OR  DRIVER  PROGRAM 


4 


CALL  DATA IN 
CALL  TAPSET 
CALL  ETGEO 


CALL  ETGDM(DTIME) 


100 

101 


200 

201 

202 

203 

204 

205 

206 

207 

208 

209 

210 
,211 

212 


WRITE  (LW, 208) 
WRITE  (LW.  209) 
WRITE  (LW. 210) 
WRITE  (LW. 211) 
WRITE  (LW. 212) 


RO.  Rl,  ROD 

(AMP(K),  DA ( K ) »  DB(K>,  DC  <K),DAD(K),  K=l,  NS) 
DTI  ME 

( VR ( I ) , V I ( I ) , DAZ ( I ) , DEL ( I ) , I = 1 . NT AP ) 

PEFF 


STOP 


FORMAT (3F 10.  1) 
FORMAT ( 15) 


FORMAT (//14H 
FORMAT ( / / 1 7H 
FORMAT (//19H 
FORMAT (//14H 
FORMAT (//14H 
FORMAT (//17H 
FORMAT ( // 1 1H 
FORMAT  (//UH 
FORMAT (//15H 
FORMAT ( //22H 
FORMAT (//11H 
FORMAT ( // 19H 
FORMAT (//l OH 
END 


XO.  YO,  ZO . /  (20X3F12.  6)  ) 

XOD.  YOD,  ZOD . /  ( 20X3F12.  6)  > 

PSI.  THETA.  PHI . /  ( 20X3F12.  6)  ) 

BP,  BQ,  BR . /  (20X3F12.  6) ) 

XM,  YM,  ZM . / (20X3F12.  6)  ) 

XMD,  YMD,  ZMD . /(20X3F12.  6)) 

NSCAT . /20XI12) 

X,  Y,  Z . /  (20X3F12.  6) ) 

RO,  Rl,  ROD . / ( 20X3F12.  6)  ) 

AMP,  DA.  DB,  DC,  DAD . /(20X5F12. 

DTI  ME . /20XF12.  6) 

VR,  VI ,  DAZ,  DEL . /  (  20X4F1 2.  6 )  ) 

PEFF . /20XE12.  5) 


6)  ) 


oftoooonnoooonoonnooonnnnonnooonnonnooonftoooonoono 
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SUBROUTINE  ETGEO 

TRANSFORMATION  TO  RADAR  SPACE  FOR  N-POINT  SCATTERER  MODEL 

IN  THIS  SUBROUTINE  WE  BEGIN  WITH  THE  MODEL  OF  AN  EXTENDED  TARGET  AND 
THE  ENGAGEMENT  GEOMETRY  IN  ORDER  TO  COMPUTE  THE  AMPLITUDE  AND  RADAR 
COORDINATES  FOR  EACH  SCATTERER  IN  THE  MODEL. 

THE  MODEL  IMPLEMENTED  IS  THE  SO-CALLED  MEDIUM-RANGE  MODEL  <SEE  MR I 
REPORT  132-44). 

ASSUMPTIONS  AND  LIMITATIONS . 

1.  ALL  SCATTERERS  ASSUMED  TO  BE  ILLUMINATED  BY  SAME  TRANSMIT 
ANTENNA  GAIN. 

2.  TARGET  ASSUMED  TO  BE  WITHIN  LINEAR  REGION  OF  MONOPULSE 
RECEIVE  BEAM. 

3.  THE  DOPPLER  SHIFT  OF  THE  TARGET  CG  IS  IMPLEMENTED  BY  MEANS 
OF  A  F I NELY-CONTRQLLABLE  DELAY  LINE  (THE  LASER  DEVICE), 
PLUS  THE  USE  OF  THE  FREQUENCY  SYNTHESIZER. 

4.  ONLY  ONE  PHYSICAL  TARGET  IS  SIMULATED  PER  CALL. 


ALL  COMMUNICATION  TO  AND  FROM  THIS  SUBROUTINE  IS  THRU  COMMON. 


ON  INPUT 


/Tl/  XO,  YO,  ZO 

XOD,  YOD,  ZOD 
PSI, THETA.  PHI 
BP.  BQ,  BR 


=  TARGET  CG  IN  INERTIAL  COORDINATES 
=  TARGET  CG  RATE  IN  INERTIAL  COORDINATES 
*  YAW,  PITCH,  ROLL  ANGLES 
■  YAW,  PITCH,  ROLL  ANGLE  BODY  RATES 


/T3/  XM.  YM,  ZM 

XMD. YMD,  ZMD 


=  MISSILE  CG  IN  INERTIAL  COORDINATES 
»  MISSILE  CG  RATE  IN  INERTIAL  COORDINATES 


/T4/  NSCAT  ■  NUMBER  OF  SCATTERERS  IN  TARGET  MODEL 
ST  »  APPROXIMATE  PHYSICAL  SIZE  OF  TARGET 
AMPMIN  *  AMPLITUDE  THRESHOLD  FOR  SCATTERERS 
X, Y, Z  *  ARRAYS  CONTAINING  SCATTERER  LOCATIONS  IN  TARGET 
COORDINATES 

/TS/  NTAP  =  NUMBER  OF  TAPS  IN  TAPPED  DELAY  LINE 
DRTAP  »  SPACING  BETWEEN  TAPS  (RANGE) 

,ON  OUTPUT . 

/T2/  CPSI »  SPSI,  .  .  .  ETC  »  SINES  AND  COSINES  OF  TARGET  ANGLES 

/T7/  NS  -  NUMBER  OF  SCATTERERS  VISIBLE 
RO  -  RANGE  TO  TARGET  CG 

Ri  -  RANGE  TO  FIRST  TAP 


oooouoooouuouoooooouooooooooouooooooooooo  u  o 


ROD  =  RANGE  RATE  OF  TARGET  CG 

AMP < J )  »  AMPLITUDE  OF  U-TH  SCATTERER 

DA(  J)  =  INCREMENTAL  A-VECTOR  OF  J-TH  SCATTERER 

DB(J>  ■  INCREMENTAL  B -VECTOR  OF  J-TH  SCATTERER 

DC(J)  =  INCREMENTAL  C-VECTQR  OF  J-TH  SCATTERER 

DAD < J >  ■  INCREMENTAL  A-VECTOR  RATE  OF  J-TH  SCATTERER 

THE  TARGET  CG  AND  MISSILE  CG  COORDINATES  ARE  IN  AN  INERTIAL  COORDINATE 
SYSTEM  REFERENCED  TO  THE  GROUND  (XY-PLANE  PARALLEL-  TO  GROUND.  2-  DOWN) 

THE  ABC-VECTORS  ARE  DEFINED  AS . 

A  -  FROM  THE  TARGET  TO  THE  MISSILE 

B  -  PARALLEL  TO  THE  GROUND,  TO  THE  LEFT  AS  VIEWED  FROM  MISSILE 
C  -  PERPENDICULAR  TO  A  AND  B  IN  RIGHT-HAND  COORDINATE  SYSTEM 

THE  TARGET  COORDINATES  ARE . 

X  -  TARGET  LONGITUDINAL  AXIS,  POSITIVE  IN  DIRECTION  OF  NOSE 
Y  -  IN  DIRECTION  OF  RIGHT  WING 
Z  -  DOWN 

THE  BODY  RATES  ARE  DEFINED  AS . 

BP  -  CW  ROTATION  RATE  ABOUT  TARGET  X-AXIS 

BQ  -  CW  ROTATION  RATE  ABOUT  TARGET  Y-AXIS 

BR  -  CW  ROTATION  RATE  ABOUT  TARGET  Z-AXIS 

THE  DIRECTION  OF  ROTATION  IS  DEFINED  LOOKING  OUT  FROM  THE  COORDINATE 
ORIGIN. 

SEE  SUBROUTINE  XFORM  FOR  A  DEFINITION  OF  THE  YAW,  PITCH,  AND  ROLL 
ANGLES. 

THE  RFSS  CHAMBER  COORDINATES  ARE  ASSUMED  TO  BE  PARALLEL  TO  THE  ABC- 
VECTORS.  RANGE  IS  IN  -A  DIRECTION,  RIGHT  AZIMUTH  IN  -B  DIRECTION,  AND 
UP  ELEVATION  IN  -C  DIRECTION. 

ALL  DISTANCES  (INCLUDING  WAVELENGTH)  MUST  BE  IN  THE  SAME  UNITS.  ALL 
ANGLES  MUST  BE  IN  RADIANS. 

< 

DIMENSION  A<3)>  B<3)>  C<3),  WA ( 3 ) ,  WB ( 3 ) ,  WC ( 3 ) 

COMMON  /Tl/  XO,  YO,  ZO,  XOD,  YOD,  ZOD,  PSI,  THETA,  PHI,  BP,  BG,  BR 

COMMON  /T2/  CPSI, SPSI, CTHETA, STHETA, CPHI, SPHI 

COMMON  /T3/  XM,  YM,  ZM,  XMD,  YMD,  ZMD 

COMMON  /T4/  NSC AT,  ST,  AMPMIN,  X  <  20 ) >  Y( 20) »  Z ( 20) 

COMMON  /T5/  NT AP, DR TAP  < 

COMMON  /T7/  NS, RO,  R 1 ,  ROD, AMP ( 20 ) , DA( 20) , DB (20 ) ,  DC ( 20 ) ,  DAD  <  20 ) 

COMPUTE  SINES  AND  COSINES  OF  ANGLES 


o  o  o  oonoooo  ooo  ooo  r»oo  non  o  o  o 
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CALL  SINCOS(PSIi  SPSI, CPSI ) 

CALL  SINCQS< THETA,  STHETA,  CTHETA) 

CALL  SINCOS(PHI,  SPHI, CPHI > 

COMPUTE  RANGE  TO  TARGET  CG  AND  A-VECTOR 

A<l>*XM-XO 

A<2)=YM-Y0 

A<3)*ZM-Z0  , 

RO-SQRT ( A ( 1 ) **2+A ( 2  >  **2+A ( 3 ) **2  > 

A(l>=A(l)/RO 

A<2)»A(2)/R0 

A(3)*A(3)/R0 

COMPUTE  RANGE  TO  FIRST  TAP 

R1-RO-.  5*<NTAP-1 )*DRTAP 

COMPUTE  RANGE  RATE  OF  TARGET  CG 

ROD»A ( 1 ) * ( XOD-XMD ) +A  <  2  >  *  <  YOD-YMD  >  +A  <  3  >  *  < ZOD-ZMD ) 

COMPUTE  B-  AND  C-VECTORS 

RHO-SQRT  <  A  <  1 ) **2+A  <  2 ) **2 ) 

B<1)»-A<2)/RH0 
B<2)»  AUJ/RHO 
B(3)=0. 

C ( 1 ) =— A  <  3 ) *B ( 2 ) 

C<2)»  A<3)*B<1> 

C<3)»RH0 

TRANSFORM  A-,  B-.  AND  C-VECTORS  TO  TARGET  COORDINATES 

CALL  /FORM < A*  WA) 

CALL  XFQRM<2, WB ) 

CALL  XFORMfC, WC> 

COMPUTE  TAROET  ASPECT  ANGLE  < ALPHA*AZ IMUTH,  BETA=ELEVATION» 

ANGL« ANCLE  TO  ROLL  AXIS) 

ALPHA-ATAN2 ( WA  <  2 ) ,  WA  < 1 ) ) 

SBETA— WA<3) 

BETA*ATAN2 < SBETA,  SORT ( 1 .  -SBETA**2 ) ) 

ANCL-ATAN2  <  SORT  < 1 .  -WA  < 1 ) **2 ) ,  WA  < 1 >  > 

LOOP  OVER  SCATTERERS 


L-l 


o  ci  n  ana 


DO  20  K=1 <  NSCAT 
SAMP*SCTAMP(K,  ANGL) 

IFtSAMP.  L£.  AMPMIN)  GO  TO  20 
AMP(L)»SAMP 

COMPUTE  INCREMENTAL  A<  S>  C  COORDINATE 

DA<L)=X(K)*WA< 1 )+Y(K)*WA<2)+Z(K)*WA<3> 

DB(L)«X(K)*WB< 1 )+Y<K>*WB(2)+Z(K)*UB<3> 

DC<L)*X(K)*WC< 1 )+Y<K)*WC(2)+Z<K)*WC<3) 

COMPUTE  INCREMENTAL  A-VECTOR  RATE  (SMALL  ANGLES  ARE  ASSUMED) 

XKD=*  Z(K)*BQ-Y(K)*BR 
YKD=*-Z ( K ) *BP+X ( K ) *BR 
ZKD=  Y(K)*BP-X(K)*BQ 
DAD  <  L ) =XKD*WA ( 1 ) +YKD*WA ( 2 ) +ZKD*WA  <  3 ) 

L*L+1 

20  CONTINUE 
NS*L-1 
RETURN 
END 
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SUBROUTINE  ETGDM(DTIME) 


GLINT  AND  DOPPLER  MODULATION  FOR  N-POINT  SCATTER  MODEL 

IN  THIS  SUBROUTINE  WE  COMPUTE  THE  GLINT  OFFSETS  AND  MODULATION  SIGNALS 
APPLIED  TO  EACK  TAP  OF  THE  TAPPED-DELAY  LINE.  IT  IS  TO  BE  CALLED 
AFTER  ETGEO  TRANSFORMS  COORDINATES  TO  RADAR  SPACE.  IT  WILL  USUALLY 
BE  CALLED  MORE  FREQUENTLY  THAN  ETGEO.  IT  IS  ALSO  THE  BEST  SUBROUTINE 
TO  PLACE  IN  THE  AP120B. 

EXCEPT  FOR  TIME,  ALL  COMMUNICATION  TO  AND  FROM  THIS  SUBROUTINE  IS  THRU 
COMMON. 

ON  INPUT . 

DTIME  -  TIME  SINCE  LAST  UPDATE  IN  TARGEO 

/T5/  NTAP  =  NUMBER  OF  TAPS  IN  TAPPED  DELAY  LINE 
DRTAP  ■  SPACING  BETWEEN  TAPS  (RANGE) 

XL  =  WAVELENGTH 

PTCDSQ  »  PRODUCT  OF  TRANSMIT  POWER,  GAIN,  AND  SQUARE  OF 
RFSS  CHAMBER  LENGTH 

/T7/  NS  -  NUMBER  OF  SCATTERERS  VISIBLE 
RO  »  RANGE  TO  TARGET  CG 

R1  -  RANGE  TO  FIRST  TAP 

AMP <  J )  =  AMPLITUDE  OF  U-TH  SCATTERER 
DA(  J)  a  INCREMENTAL  A-VECTOR  OF  J-TH  SCATTERER 
DB(J)  *  INCREMENTAL  B -VECTOR  OF  J-TH  SCATTERER 
DAD (J)  =  INCREMENTAL  A-VECTOR  RATE  OF  J-TH  SCATTERER 

ON  OUTPUT . 

/T8/  PEFF  «  EFFECTIVE  RADIATED  POWER  AT  RFSS  ARRAY 

VR ( I )  «  IN-PHASE  MODULATION  SIGNAL  TO  I-TH  TAP 

VI ( I )  «  QUADRATURE  MODULATION  SIGNAL  TO  I-TH  TAP 

DAZ < 1 )  a  CLINT  OFFSET  (AZIMUTH)  FOR  I-TH  TAP 
DEL< I )  a  CLINT  OFFSET  (ELEVATION)  FOR  I-TH  TAP 

THE  PARAMETER  PMIN  IS  JUST  SOME  SMALL  NUMBER  TO  PREVENT  DIVIDE  BY  ZERO 

ARRAYS  VBR, VBI.  VCR, VCI  MUST  BE  DIMENSIONED  AS  LARGE  AS  NTAP 

DIMENSION  V3R<8).  VBI(S), VCR(S),  VC 1(8) 

DIMENSION  SS(4) , CC (4 ) 
t  DIMENSION  TW (4 ) 

COMMON  /T5/  NTAP, DRTAP. DRGATE,  XL,  PTGDSQ 

COMMON  /T 7/  NS,  RO, Rl, ROD, AMP ( 20) ,  DA( 20) ,  DB ( 20 ) ,  DC ( 20) ,  DAD (20) 
COMMON  /TS/  PEFF,  VR(8), VI ( 8) ,  DAZ < 8) ,  DEL (8) 

DATA  PMIN/1.  E-10/ 
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DATA  FOUR? 1/12.  5663706/ 

ZERO  ARRAYS 

CALL  XMIT<-NTAP» 0.  »VR) 

CALL  XMIT<-NTAP. 0.  , VI ) 

CALL  XMIT<-NTAP.  0.  ,  VBR) 

CALL  XMIT<-NTAP» 0.  » VBI ) 

CALL  XMIT<-NTAP.  0.  .  VCR) 

CALL  XMIT<-NTAP,  0.  ,  VCI) 

CALL  XMIT (-NTAP.  0.  *  DAZ ) 

CALL  XMIT < -NTAP. 0.  # DEL) 

LOOP  OVER  NS  SCATTERERS 

DO  40  J=l»  NS 

COMPUTE  TAP  NUMBER  OF  FIRST  TAP  (ITAP)  AND  FRACTION  (P) 

R=RO—  <  DA ( J ) +DAD  <  J ) *DT IME ) 

P*  <  R-R 1 ) /DRTAP+1 00. 

ITAP=P 
P-P-ITAP 
ITAP=ITAP— 100 

COMPUTE  RANGE  DIFFERENCE  FROM  TAP  NUMBER  ITAP 
DR»(P+1.  )*DRTAP 
FIND  TAP  WEIGHTS 

CALL  TAPWTStP.  TW) 

COMPUTE  PHASE  ON  FOUR  TAPS 
DO  20  1*1/4 

CALL  SINCOS<-FOURPI*DR/XL.  S.  C) 

SS(I)=S*AMP<J)*TW< I) 

CCC I )*C*AMP< J)*TW< I ) 

DR»DR-DRTAP 
20  CONTINUE 

LOOP  OVER  UP  TO  FOUR  TAPS  AND  INCREMENT  ARRAYS 

IF<  ITAP.  GT.  NTAP)  CO  TO  40 
IF<  ITAP.  LT.  -2)  GO  TO  40 
I 1*MAX0< ITAP# 1) 

I2-MIN0< ITAP+3,  NTAP) 

II-I1-ITAP 
DO  30  1*11.  12 
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30 

40 


II-II+l 

VR  <I)=VR  (I>+CC(II) 

VI  (I)=VI  <I)+SS(II) 

VBR ( I ) =VBR ( I ) +CC (II) *DB ( J ) 
VB 1  ( I )  =VB  Id)  +SS  (II)  *DB  ( J ) 
VCR  < I )=VCR ( I )+CC ( 1 1 >*DC ( J) 
VC  I  ( I )  =VC  Id)  +SS  (II)  *DC  ( J  > 
CONTINUE 
CONTINUE 


COMPUTE  GLINT 


OFFSETS  FOR  EACH  TAP  AND  PEAK  POWER 


PEAK=0. 

DO  50  1=1/  NTAP 
POW=VR ( I ) **2+VI < I >  **2 
IFCPOW.  GT.  PEAK)  PEAK=POW 
IF(POW.  LT.  PMIN)  GO  TO  50 

DAZ  < I ) =-<  VBR ( I ) *VR  < I ) +VB I(I)*VI(I))/( RO*POW ) 
DEL ( I ) — ( VCR ( I ) *VR ( I ) +VC I(I)*VI(I>)/( RO*POW ) 
50  CONTINUE  \ 


NORMALIZE  AMPLITUDE 


ANORM=SGRT (PEAK ) 
DO  60  1=1, NTAP 
VR  < I ) =VR ( I ) / ANORM 
VI < I )=VI < I ) /ANORM 
60  CONTINUE 


COMPUTE  EFFECTIVE  RF  POWER 

PEFF=PEAK*PTGDSQ/ (F0URPI*R0**4> 

RETURN 

END 


/ 
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SUBROUTINE  XFORM<A.  W> 

IN  THIS  SUBROUTINE  WE  TRANSFORM  A  VECTOR  (A)  IN  INERTIAL  COORDINATES 
TO  A  VECTOR  <W)  IN  TARGET  COORDINATES.  THE  COORDINATE  ROTATIONS.  IN 
THE  ORDER  OF  APPLICATION.  ARE . 

PSI  -  CW  ROTATION  OF  Z-AXIS 
THETA  -  CW  ROTATION  OF  Y-AXIS 
PHI  =  CW  ROTATION  OF  X-AXIS 

THE  DIRECTION  OF  ROTATION  IS  DEFINED  LOOKING  OUT  FROM  THE-  COORDINATE 
ORIGIN.  IN  THIS  SUBROUTINE  THE  SINES  AND  COSINES  OF  THE  ANCLES  ARE 
INPUT  THROUGH  COMMON  /T2/. 

DIMENSION  A(3).W(3) 

COMMON  /T2/  CPSI.  SPSI.  CTHETA.  STHETA.  CPHI.  SPHI 
UX=  A<1)*CPSI+A<2)*SPSI 
UY*-A ( 1 ) *SPS I +A < 2 ) *CPS I 
UZ-  A< 3) 

VX»  UX *CTHET A-UZ*STHET A 
VY-  UY 

VZ-  UX*STHETA+UZ*C THETA 

wm«  vx 

W<2>=  VY*CPHI+VZ*SPHI 

W(3)=-VY*SPHI+VZ*CPHI  , 

RETURN 

END 
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SUBROUTINE  TAPWTS(P,  TW) 


IN  THIS  SUBROUTINE  FOUR  TAP  WEIGHTS  ARE  RETURNED  IN  ARRAY  TW  ACCORDING 
TO  THE  FRACTION  P.  THE  WEIGHTS  ARE  EXTRACTED  FROM  A  PRECOMPUTED  TABLE 
(SEE  SUBROUTINE  TAPSET). 

ARRAY  TWARAY  IS  USED  AS  IF  IT  WERE  DIMENSIONED  (4, NARAY). 


DIMENSION  TW(4) 

COMMON  /T6/  NARAY.  TWARAY (1) 

DATA  LW/6/ 

INDEX“< NARAY- 1 ) *P+1.  5 

CALL  XMIT(4f TWARAY(4*IND£X-3).  TW) 


RETURN 

END 


/ 


A 
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SUBROUTINE  TAPSET 


IN  THIS  SUBROUTINE  THE  TAP  WEIGHT  TABLE  IS  COMPUTED.  IT  IS  A 
COMPANION  SUBROUTINE  TO  TAPWTS,  AND  IT  IS  TO  BE  CALLED  AS  AN  INITIAL 
IZATION  STEP  PRIOR  TO  THE  BEGINNING  OF  THE  SIMULATED  MISSION. 

/T5/  DRTAP  *  SPACING  BETWEEN  TAPS  (RANGE) 

DRCATE  *  SPACING  BETWEEN  RECEIVER  GATES  (RANGE) 

ARRAY  TWARAY  MUST  BE  DIMENSIONED  AS  LARGE  AS  4*NARAY. 

DIMENSION  A<4, 4)* X<4) 

COMMON  /T5/  NTAP. DRTAP, DR GATE 
COMMON  HU  NARAY,  TWARAY  (I ) 

D-DRTAP/DRGATE 

L*1 

DO  30  K=»l»  NARAY 
P»(K-1)/FL0AT<NARAY-1) 

DO  10  J»l,4 
X<J)*CHI<D*<P+2-J> ) 

10  CONTINUE 
DO  20  1  =  1,4 
DO  20  J=l,  4 
A<I,  J)=CHI<D*(I-J) ) 

20  CONTINUE 

CALL  SIMG< A, X, 4, I ERR) 

IF<  IERR.  GT.  0)  STOP 
CALL  XMIT ( 4,  X,  TWARAY <L) ) 

L*L+4 

30  CONTINUE 
RETURN 
END 
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FUNCTION  CHI (P) 

RANGE  GATE  RESPONSE.  THE  ARGUMENT  P  IS  THE  RANGE  MISMATCH  NORMALIZED 
TO  THE  RECEIVER  GATE  SPACING.  INTERPOLATION  IS  USED  ON  THE  SAMPLES 
STORED  IN  THE  A-ARRAY,  WHERE  THE  SPACING  IS  0.  1  UNIT. 

THE  RESIDUAL  ERROR  IN  THE  INTERPOLATION  IS  LESS  THAN  .  0003 

P  MUST  BE  LESS  THAN  1.  5  IN  MAGNITUDE. 

THE  SAMPLES  ARE  OF  THE  RESPONSE  DERIVED  IN  MR I  REPORT  149-4. 


DIMENSION  A(  18) 

DATA  A/1.  00000*  .98104,  .92193,  .81903,  .67431,  .50112,  .32385, 

1  .  17071,  .  06308,  .  00731,  -.  00651,  .  00182,  .  01262,  .  01458, 

2  .  00713.  00313,  -.  00898,  -.  00762  / 

H-10.  *ABS(P) 

IF<H.  OT.  15.  )  STOP  55 

I-H 

H=H-I 

IP1-I+1 

IP2-I+2 

IP3-I+3 

IF(  I.  LE.  0)  1-2 

CHI—.  1 66667 *H* < H—  1 .  >*<H-2.  )*A(I)+.  5*<H**2-1.  )*<H-2.  )*A(IP1> 

1  -.  5*H*<H+1.  ) * ( H-2.  )*AUP2)+.  166667*H*(H**2~1.  )*A<IP3) 

RETURN 

END 


/ 
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SUBROUTINE  SXMQ(A»  B» N» IERR) 

SOLVES  SET  OF  N  SIMULTANEOUS  EQUATIONS . 

A  *  X  -  B  SUM  ( A< I, J)*X ( J) )  »  B ( I ) 

WHERE  ARRAY  A  IS  2-DIMENSIONAL.  ARRAY  X  IS  RETURNED  IN  ARRAY  B.  AND 
ARRAY  A  IS  DESTROYED.  COMPUTATION  IS  VALID. IF  IERR=0 

DIMENSION  A(  1 ).  B(  1 ) 

IERR  *  0 

IF  (N.  ©T.  O)  ©0  TO  10 

IERR  «  1 

RETURN 

FORWARD  SOLUTION 

10  TOL  ■  0.  0 
KS  «  0 
JJ  *  -N 
DO  65  J  »  1/  N 
JY  «  J  +  1 
JJ  -  JJ  +  N  +  1 
BICA  -  0. 

IT  »  JJ  -  J 
DO  30  I  *  J.N 

SEARCH  FOR  MAXIMUM  COEFFICIENT  IN  COLUMN 
IJ  »  IT  +  I 

IF  (ABS(BIGA)— ABS< A( IJ) > )  20.  30.  30 

20  BICA  »  A( I J) 

IMAX  ■  I 
30  CONTINUE 

TEST  FOR  PIVOT  LESS  THAN  TOLERANCE  (SINGULAR  MATRIX) 

IF  <ABS< BICA) -TOL)  35.  35,  40 
35  IERR  -  2 
RETURN 

INTERCHANGE  ROWS  IF  NECESSARY 

40  II  -  J  ♦  N* < J— 2 ) 

IT  -  IMAX  -  J 
DO  50  K  -  J,N 

11  *  II  +  N 

12  •  II  ♦  IT 
SAVE  -  A< I 1 ) 

A( II )  •  A< 12) 
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r  ; 


»  -  A... 


A( 12)  *  SAVE 

DIVIDE  EQUATION  BY  LEADING  COEFFICIENT 

50  A( II )  »  A< I 1 ) /BIGA 
SAVE  =  B<IMAX) 

B(IMAX)  =  B<J) 

B(J)  *  SAVE/BIGA 

ELIMINATE  NEXT  VARIABLE 

IF  <J-N)  55,  70,  55 

55  IQS  »  N*<J-1> 

DO  65  IX  -  JY, N 
IXJ  =  IQS  +  IX 
IT  -  J  -  IX 
DO  60  JX  *  JY,  N 
IXJX  *  N*< JX— 1 )  +  IX 
JJX  «  IXJX  +  IT 

60  A(IXJX)  »  A(IXJX)  -  <  A  < I X J ) *A  <  JJX  > ) 

65  B(IX)  «  B < IX )  -  <B<J)*A<IXJ>> 

BA  SOLUTION 

70  NY  -  N  -  1 
IT  «  N*N 
DO  80  J  -  1, NY 
IA  -  IT  -  J 
IB  -  N  -  J 
IC  «  N 

DO  80  K  -  l.J 

B( IB)  «  B < IB )  -  A( I A) *B ( IC ) 

IA  a  IA  -  N 
80  IC  a  IC  -  1 
RETURN 
END 
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SUBROUTINE  XMIT  <N,  A,  B ) 

« 

IN  THIS  SUBROUTINE  WE  EITHER  TRANSMIT  ARRAY  A  TO  ARRAY  B  (IF  N.  GT.  0> 

OR  WE  TRANSMIT  THE  CONSTANT  A  TO  ARRAY  B  (IF  N.  LT.  0).  IN  EITHER  CASE 
THE  ARRAY  LENGTH  IS  IABS(N>. 

j 

THIS  SUBROUTINE  SHOULD  BE  WRITTEN  IN  ASSEMBLY  LANGUAGE  ^ 

DIMENSION  A(  1 ) ,  B ( 1 ) 

IF(N)  10,, 20,  25 
10  NN=-N 
AA-A(l) 

DO  15  K“l» NN 

B(K)=AA  i 

15  CONTINUE 
20  RETURN 
25  DO  30  K=l,  N 
B(K)»A(K) 

30  CONTINUE 

RETURN  ^ 

END  J 


/ 


a 
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#  SUBROUTINE  SINCOS< ARG. S,  C> 

THIS  SUBROUTINE  SHOULD  BE  WRITTEN  IN  ASSEMBLY  LANGUAGE,  USING  THE 
TABLE-LOOKUP  METHOD  DESCRIBED  BY  MITCHELL  (RADAR  SIGNAL  SIMULATION). 

S-SIN(ARG) 

•  C-COS(ARG) 

RETURN 

END 
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SUBROUTINE  ETGD1 (DTIME) 

GLINT  AND  DOPPLER  MODULATION  FOR  N-POINT  SCATTER  MODEL 
NO  RANGE  EXTENSION 
SUBROUTINE  REPLACES  ETGDM 
ON  INPUT . 

DTIME  -  TIME  SINCE  LAST  UPDATE  IN  TAR GEO 
/T5/  XL  -  WAVELENGTH 

PTGDSQ  »  PRODUCT  OF  TRANSMIT  POWER,  GAIN.  AND  SQUARE  OF 
RFSS  CHAMBER  LENGTH 

/T7/  NS  *  NUMBER  OF  SCATTERERS  VISIBLE 

AMP(J)  =  AMPLITUDE  OF  J-TH  SCATTER ER 
DA( J)  »  INCREMENTAL  A-VECTOR  OF  J-TH  SCATTERER 

DB(J)  «  INCREMENTAL  B -VECTOR  OF  J-TH  SCATTERER 

DC(  J)  -  INCREMENTAL  C-VECTDR  OF  J-TH  SCATTERER 

DAD ( J )  ■  INCREMENTAL  A-VECTOR  RATE  OF  J-TH  SCATTERER 

ON  OUTPUT . 

/T8/  PEFF  *  EFFECTIVE  RADIATED  POWER  AT  RFSS  ARRAY 

/T9/  VR.  VI  -  DOPPLER  MODULATION  SIGNAL 

DR.  DAZ. DEL  -  RANGE, AZIMUTH,  AND  ELEVATION  GLINT  OFFSETS 

COMMON  /T5/  NTAP. DRTAP.  DRGATE,  XL,  PTGDSQ 

COMMON  /T7/  NS,  RO, Rl, ROD, AMP  <20) «  DA( 20) ,  DB ( 20 ) » DC (20) ,  DAD (20) 

COMMON  /T8/  PEFF 

COMMON  /T9/  VR, VI, DR, DAZ, DEL 

DATA  FOURPI/12.  5663706/ 

ZERO  ACCUMULATORS 

VR-O. 

VI -O. 

VAR-O. 

VAI— O. 

VBR-O. 

VBI-O. 

VCR-O. 

VCI-O. 

LOOP  OVER  NS  SCATTERERS 


DO  40  J-l.NS 
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#  CALL  SINCOS<FOURPI*<DA< J)+DAD< J)*DTIME) /XL#  S»  C) 
C»C*AMP<J) 

S*S*AMP<J> 

VR  *VR  +C 
VI  *VI  +S 
VAR»VAR+C*DA<J) 

•  VAI*VAI+S*DA<J) 

VBR=VBR+C*DB<J) 

VBI=VBI+S*DB<J) 

VCR*VCR+C*DC<J> 

VCI*VCI+S*DC<J) 

40  CONTINUE 

P0W=VR**2+VI**2 
AMPL»SGRT<POW> 

COMPUTE  CLINT  OFFSETS 

DR  — <VAR*VR+VAI*VI)/POW 
DAZ»-<VBR*VR+VBI*VI)/<RO*POW) 

DEL*- <  VCR  *VR  +VC I *V I ) / <  RO*POW ) 

COMPUTE  EFFECTIVE  RF  POWER 

PEFF»POW*PTGDSQ/ <  FOURP I *R0**4 ) 

NORMALIZE 

VR*VR/AMPL 
VI-VI/AMPL 

RETURN 
END 
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SUBROUTINE  DATAIN 

READS  TARGET  SCATTERING  DATA  SUPPLIED  BY  M.  MUMFORD  (SEE  SCTAMP) 

DIMENSION  1 A  ( 1 ) .  A A ( 4 ) »  X  X  ( 4 ) »  Y Y ( 4 ) .  22(4) 

COMMON  /DP/  P(IOO). IP<100> 

COMMON  /DG/  0(918) 

COMMON  /T4/  NSCAT 
DATA  LR, LW/5. 6/ 

NSCAT* 10 
M-l 

DO  20  1*1 4 NSCAT 
PRINT  99.  I 
L«10*(I-1> 

10  L-L+l 

READ  <LR.  100)  IA.  P  (L) .  AA,  XX.  YY.  22 
WRITE  (LW,  100)  IA.  P(L).  AA.  XX.  YY.  22 
IP (L)*M 
IA-IA-2 

CALL  XMIT<17. IA.Q(M)) 

M-M+17 

IF(P (L) .  LT.  ISO.  )  GO  TO  10 
20  CONTINUE 
RETURN 

99  FORMAT </29H  TARGET  DATA  FOR  SCATTERER  N0I3//) 

100  FORMAT (  1X1 1.  12XF8.  3.  4E14.  8/(22X4E14.  8)  ) 

END 


/ 
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FUNCTION  SCT AMP ( K,  ANGL ) 

v 

IN  THIS  SUBROUTINE  WE  COMPUTE  THE  AMPLITUDE  (SQRT(RCS))  OF  THE  K-TH 
SCATTERER  AS  VIEWED  FROM  THE  TARGET  ASPECT . 

ANGL  =  ANGLE  FROM  ROLL  AXIS  MEASURED  FROM  NOSE  (RAD) 

IN  ADDITION  IN  COMMON  /T4/ . * 

ST  *  BIAS  THAT  IS  ADDED  TO  ANGL  (RAD) 

THIS  SUBROUTINE  ACCESSES  TARGET  DATA  SUPPLIED  BY  MIKE  MUMFORD  AT  NWC/ 
CHINA  LAKE  IN  THE  FORMAT  DEFINED  BY  A  COMPUTER  PROGRAM  WRITTEN  5/11/78 
BY  E.  HUTTON  X3219. 

DIMENSION  IA(1). AA<4), XX (4), YY<4), ZZ(4) 

COMMON  /T4/  NSC AT »  ST  *  AMPMIN»  X  <  20 ) ,  Y  <  20) ,  2 ( 20) 

COMMON  /DP/  P( 100), IP(IOO) 

COMMON  /DQ/  Q<918) 

ANG=ABS ( ANGL+ST ) *57.  2957795 
IF  (ANG.  GT.  180.  )  ANG*180. 

I1*10*(K-1)+1 
12*11+9 
DO  20  1*11,  12 
IF  (ANG.  LT.  P  <  I )  )  GO  TO  25 
20  CONTINUE 
25  M-IP(I) 

CALL  XMIT (17,Q(M),  IA) 

I'(IA)  30,31,32 

30  SCTAMP=AA ( 1 ) +ANG* < AA ( 2 ) +ANG* ( AA ( 3 ) +ANG*AA ( 4 ) ) ) 

GO  TO  35 

31  SCTAMP=EXP(AA( 1 )+ANG*AA(2) ) 

GO  TO  35 

32  SCTAMP =EXP ( AA ( 1 ) +ANG*AA ( 2 ) ) -EXP  <  AA  <  3 ) +ANG*AA ( 4 ) ) 

33  SCTAMP*.  09004*SCTAMP 
IFCSCTAMP.  LT.  AMPMIN)  RETURN 

X (K)«XX ( 1 ) +ANG*( XX ( 2 ) +ANG*( XX ( 3 ) +ANG*XX (4 ) ) ) 

Y ( K ) *YY ( 1 ) +ANG* ( YY ( 2 ) +ANG*  <  YY ( 3 ) + ANG*YY  <  4 ) ) ) 

Z  <K)=ZZ ( 1 )+ANG*( ZZ (2)+ANG*<  ZZ  <  3)+ANG*ZZ (4) ) ) 

Y(K)— Y<K) 

Z<K>—  Z<K) 

RETU8N 
'  END 


